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ABSTRACT 


The increase in trans-oceanic travel and maritime operations have 
created a need for new and improved search and rescue techniques and 
facilities. This study briefly describes the various phases of search 
and rescue operations and reviews the doctrine currently used by this 
Nation as set forth in the National Search and Rescue Manual, 

Attention is then directed to developing a computer program which 
utilizes current doctrine in the solution of one portion of the search and 
rescue problem. The program determines the most probable location of 
survivors (datum) and the associated search areas at any specified 
time for an off-shore distress case with an estimated reported distress 
incident position. The resulting program, written in Fortran 60, is 
. included with sample test data. 

Finally, recent developments in computer-generated informations 
systems which may offer improvements in the existing doctrine are 


discussed, 


il 








TABLE OF CONTENTS 

CHAPTER PAGE 
bs STATEMENT OF THE PROBLEM) .0.0 stecienesneuatersacr-icuencnenenemer erm 
Theoretical Background......... oS Spelene o eceenemene ore tepeus ee 
Statementios Special Purposes < o:c:s.s:< cca) ese scorer Chereiaetele nan 
TN CBEIOI MSI. 5 «0 os 0 0 6.46 « © srotelibesers 6 ¢ a lebettienenelebelsl esstensienenenmmnE) 
Definitiomssor Key Terms... sso» scmss olen cs 66 co cleteericn nee 
Times SOIC Cis iclsislisivieis ¢ « o bic swe 0.0 « «614 0 oils ft ticle eleietermanne” 
DIGGERS SS COLUM. ic siete isl cele © eteroLbhi te cl bicker ohobelarsie « ole) aiercremma) 
Reseuc  Gocrdination Genter. . . «. stele s samme « « clmeieneiet anne 
Oi CemeRCOMMANGCT. ¢-c. « «+ 6 es so 000 siele cis ee cle siete sine 
Osc Whi st eer eyo SoCo cou | 
Distresseincidentspositlon. < « sis.cis0%.s <. s.sieieies wo cts Peer? 
BViuQerirre re te mes hs Op G. 
DCARGM ATCA. 4s oe ‘shyre 4 4:0 4:0. 0:0 olemehille c/s 10s 0 oan SUonSRenaNeNeneT ORC mmm 
SCArChemaGUs’s . Js oo 0 Wout roc © eucte gh cietciere s/cleiseetey enone MERE 
Total probable error of pOSitiOoNn. «sss ties cose cc ee Bad 9 
Ii Ee POSITION CTVOK, ace o:o:0ce:einecerancieneretebenels vicle-« sichetele on tat 
Navigation error of the search unit..cccccccccvccces 9 
Life *vabtwelyert CUPOL tes) ws omeeielel cleo ele herein Cielele + oldies +o almene 
Drifitcogs ax S's cc 0 e's 0 + otece ef. eerqrcvetebanede: 4. otetere.ci cis ove. etenenencrn Dane 
Drift (tage) W sco os so + 0 5 s Gipierers aiemetenenerc a aieielere sane) 
Drdit Ummm). ps0 ses 66 oldie apie ocavoustutore s/o oy Sier stole cathe Deen! 


iii 








CHAPTER PAGE 

Local wind CUrrenty 5 ciciors oo chcsoneleisiolel sis intsleusiciy oheraicl creme 

LG @WaliVictweienaue.s si bie isic © oisi sities ol bis tenehe be oh bieieie ote sist terran Cooke 

Parachute drift. ..a <5 oss sible © @ Sache «elec cleissoret bin) eie ter tieec ren emmee 
ASSUMPPLONSS « ou 0c: Susice5 oS opepelene, wwe aikereisireic ts -usiene enone enn eS 
TEUIMUITALLONS s. sicisve 0.0 6 Sierete se se) o)o etensl vilcsonctie lobeNelcl ether ienenene nema: 
Research Significance suseswc s 6 -arse 0 sl ctekelcs 0.410 e)encrener owe ererer nD 

r. REVIEW O@F THE LITERATURE sc co. 3:0 < scclpeletar pice <cclep bol tee enennles 


Doctrine on the Determination of the Search Area 
According to the National Search and Rescue Manual... 18 


Studies Dealing with Parachute Drift and Ocean 
Current WL teekéso Soe @ 6) 0ce 0 0 eo te OIeIS ee 32 


Parachute rift... & s:e-s:0 «oe veueuegaiawenel baeaeueuens ¥ cecee eieicsars Ge 
Ocean current OTE og nce ok ele eee 35 


Note on Additional Literature Relating to the Overall 
Search and”Rescue Problem...seeceescboceccsssecttet JG 


ele DELEON areiacetaite sieve. 9) oie alot) eulerotetenetelenetetsl state ste enoiemenel too 7 men 
IWVS TIO Clip saya kiotsie ist Bice ssi ololelencltel cle eleucteteleiele <icuelotenerstsieicn sist mmmnT 
FPLrOCCAUIC. Fos cclce ss bisle @ eunaletonesiot EU NEI. © 6G Biel c pisisaie ee 
Determination of datum at any Specified time........ 42 

ME CWAY «os s:0i tiene clone oeiene § one cueneiere of tc ciel el el tie sieieteno yn aa 

Local wind CUment Sora a steele alates a c's cls 6+ as «oleae nee 

Average SC€Q CUITENTIS . . occ chit h 105000 0 ee ocle ser tu 

Resultant total: drilte, . <tutrpeesstereisst » occ s 6 cleterctanerenene me 


iv 








CHAPTER PAGE 


Determination of total probable error of position.....0. v4 
Initial position error of the distress incident......00.5 J9 


Navigational error of the search unit. ...c0c0c00%8000e0e0 Vo 


coy, ne so ee er ee - 61 


Determination of the limits of the search area...o.evcoe Gl 


niV, RES U Lib Grew sc eG ccc ete eo cc 6.05 ble 6 se ememenel ome eee 20 


eat. 


Taal Dats @CoOmpusatiOmes. oc... vos 60860 0.0 0 0 GEKe oo cme ame 


V, SUMMARY AND CONCLUSIONS. .ccccscccscecctccstevscs JB 


Summary, o@uvu 9300 vodgRs @OOoOsodosedaoesdreddgedoddscgs Ggeveoeetoev Goegsoev 9 76 
Conclusions ooveoceocoecesectcgcseaceecetetaoegoedooeoaoesecoet oaeeaoee 53 639 8 7 6 
49 


Pr GUO NS oi 5.55 25 oars, 65 2G Sree eat Breaks. © 20 et Oe ee 


Recommendations of Areas for Further Study.....0,-0%e00e0¢+ 80 


Determination of a datum and associated search 
area wken the time of the incident and estimated 


PESiHOn aTe NOt KNOWN: . . occ 2-c coc 36 ciwise o cleric a 


Development of a fully operational free fall and 


parachute drift PFOGTAM. 500040072000 09000000000900000 80 


Development of a program to handle small craft 
tienore) dictres@ CASES. ...c0cco rs esccccocecs sco Mpa 


Development of a program relating search area, 
s@arch units, and search patterns... ccsccccorccccce OL 


studies evaluating the frequency of distress cases 
and the location of SAR facilities......c0cceccccccee O2 


C@ntivitued drift analysi&...cccocccech oo ckcco esr osc ces MOM 


BRB LIOGRAPHY occs. cc 6 ou eae oe ee awe a Ree cca cco eemE Ce eee teen ene ee aeeeeernee 84 


APPENDIX A. Flowchart Symbols and DefinitionsS...cocsccccscees 89 


VW 
st S 3 ey is es ’ v © 9 


oar & 








APPENDIX B. 


APPENDIX CG. 


APPENDIX D. 


APPENDIX E. 


APPENDIX F, 


Program SAR oeee7e37eoeee#e¢@ eoeeoeecdceee?ee@@ 8 @ @ @—UOhUSHUOCUO oee¢8 8 @ 91 
"Parachute 1D) gh 4 hae ne 8 Re ee a ae Wes 


Extracts from the Study: "Free Fall and 
Parachute MAE «6 6 a hie ackieeie-0 aC on eee 156 


“Computer Produced Synoptic Analyses of 
Surface Currents and Their Application for 


INEM SEO easiest os 4; 6 0 00 sa% 4 0-0 « Ciswes 6 ew a. erenacLonenameeeene 174 


Suriace Drift Ellipse Program. ....e.sceccccaceccse «clean 


Vi 





LIST OF PIGURES 
. FIGURE PAGE 


lee Components of Drift, D, When Flucmations. are Small 
or Existing Conditions are known with a High Degree 
Of Agigdemne ts @ cgi’ hke comes. 6 Semmes oor aie 0 Ss ames 6.5626 0.6 eT 


. ee A Minimax Plot of Drift Components when Fluctuations 
are Large or Uncertainty Exists as to the Exact 
Conditions ar the MERI 0 0 SOOT DORE DR USER SEIS cou 20: i] 


ae ae liga. Ses, a. mee 


3, Cracniedl Solution of Drift sacl verre). sie oo 0 eles elo eae 


4. Repeated Expansion Search Areas for Stationary Datum 

PO aniline ete hrn: a SoSca es oo eared Reais eva easton BM AVE 6 cBsero oe Oar eener ee 3h 
De Repeated Expansion Search Areas for Moving Datum Point... 31 
6. The Basic Process of Determining the Search Area..cccerce 41 


,, %7A. The Factors Associated with the Forecasted Datum......6. 43 
7-B. Flow Diagram for Determining Leeway Drift Components... 47 


_ 7@-C. Flow Diagram for Determining Local Wind Current 
DPT Demi oe Cite vc 4 0 + siete 0 60000 dog ose scoc gens s MES es ° 50 


« 7°D, Flow Diagram for Determining the Average Sea Current 
Dye fits GOWROMMITE., « Qatsie cc ocdeens clehe co ude 0c cc 00 + o ctutelems SL 


: WeR-E. Flow Diagram Combining Drift Vectors and Datum, 
“ Yielding Forecasted Datumtjscccesceeccecccccccscacce ol 


' 8. Complete Drift Picture for Life Raft. @eogoogooev 0o0v; @oo0o0gcedodegedgedgdvce@#e god & 52 
§ ° Complete Drift Pieture for Life Boat. @eooogoeve@egcdgeeeeeeese oe Oe oD 92 
a 10-A, Components of Total Probable Error of Position.....c..csees OF 


19-B. Flow Diagram for the Determination of the Initial 
Position Error, arr ee ber re geen ren = 4 fe 57 


10-C,. Flow Diagram for the Determination of the Navigational 
Error of the Search Unit, NY, 6 bcei@e eS 6 ee boone a oe tone eee od 


16-D, Flow Diagram for the Determination of Drift Error, de...... 62 


Vis 
. t 
pets oy pte est d& . d. , “ ‘ gf I as 1 ok 








FIGURE PAGE 


10-E. Flow Diagram for the Final Determination of Total 
Probable Error of Position, Ce iw es 6 6 © clele ele shapale eo elclete lela sete 63 


Ike Flow Diagram for the Determination of the Search Area..... 64 


< 


2. General Flow Diagram for Equivalent Temperatures at 
Various Wind VEIOGCIEICS . . < 6 sccétemeetes o SE ee 67 


nS. Flow Diagram for Immersion Hypothermia......ccccccoecee 09 


vill 








CHAPTER I 


oTATEMENT OF THE PROBLEM 


Man's continuing struggle against the elements at sea is noted 
. throughout history. In the Book of Psalms it is recorded: 


They that go down to the sea in ships, that do business 
in great waters; 

These see the works of the LORD, and his wonders in 
the deep. 

For he commandeth, and raiseth the storm wind, which 
lifteth up the waves thereaf. 

They meunt up to the heaven, they go down again to 
the deptns: their soul is melted because of trouble. 

They reel to and fro, and stagger like a drunken man, 
and are at their wit's end. | 


Each wear literally thousands of peop}e find themselves in distress 
While located upon the high seas or upon waters over which the United 
States has jurisdiction. Fortunately for most of them, their rescue is 
only a matter of waiting a relatively short time. This life-saving 
assistance usually appears in the form of one or more units operated by 
the United States Coast Guard. Specifically charged by Title 14, 
United States Code, * the United States Coast Guard carries the 


responsibility of establishing, operating and maintaining these facilities 


oa 


The Holy Bible (King James Version; Cleveland: The World Publishing 
Company, n.d.), Psalms 107:23-27. 


2United States Code, Title 14 (Coast Guard Chapter 1, Section 2, 
Vol, III, (1958 edition); Washington: Government Printing Office, 1959), 


p. 2286. 


a 








necessary to provide such assistance as may be required by any person 


or persons, vessels or aircraft in distress upon these waters. 


I, THEORETICAL BACKGROUND 


The task of successfully conducting search and rescue operations 
at sea has been a challenge ever since man first ventured out upon the 
hydrosphere. Modern technology has contributed greatly to providing 
- the answer to this challenge. Nevertheless, each individual case 
presents the would=be=rescuers with many complex problems faught with 
uncertainties and the capriciousness of Nature. 

To increase man's effectiveness in rendering assistance to his 
fellow man on the high seas, the martime nations of the world have 
combined their efforts in developing and improving search and rescue 
techniques, establishing necessary Communication systems, and 
providing for mutual assistance and cooperation in large scale operations. 
The National Search and Rescue Manual, CG-308, developed by the 
United States Coast Guard, serves as this Nation's leading guide in 
this area, 3 


3United States Coast Guard, National Search and Rescue Manual, 
CG=308 (Washington: Government Printing Office, 1959, Amended 1964). 








Generally, search and rescue operations can be divided into four main 
phases. The first of these is the "Alert" when a distress or potential 
distress situation is suspected. This would occur, for example, when 
an aircraft or boat is overdue or has failed to meet routine scheduled 
communications checks. During this phase efforts are directed to 
determine if, in fact, there is an overdue craft, and whether or not 
communications can be established with the craft in question. 

The second, the "Distress" phase, follows when the "Alert" fails 
to establish the fact that a distress does not exist. Here one finds 
efforts being directed to determine the most likely location of the 
distressed unit with an associated area for searching and the assignment 
and dispatching of search units to the scene of the search. 

The third is the "Search" phase wherein those search units previously 
designated conduct searches of the area determined in the second phase. 

The fourth and last phase is the "Rescue" which hopefully follows 
a successful "Search." 

It is evident when a distress incident does exist that no amount of 
searching will lead to rescue if this searching has been conducted in an 
area which does not contain the unit in distress. Therefore, it is 
essential that the second phase produce sound information and directions 


for the search units. This requires the careful analysis and evaluation of 








many complex and interrelated factors. Some of these factors are the 
nature of the distress, reliability of navigational information available 
on.the probable position of the distressed unit, the effects of ecean currents, 
local wind-driven currents, sea and swell conditions, life raft drift or 
leeway, possible effect of parachute drift if such is present in the 
distress incident under consideration, estimated time of survival under 
conditions existing at the scene, time involved for the search units to 
reach the scene with the associated change in site location caused by 
the factors stated earlier, the navigational errors of the search units, 
and the various probabilities associated with detection in different size 
areas. 

At the present time the evaluation and analysis of these factors are 
Carried out by human decision-makers with the aif! of excellent guide- 
lines, charts, graphs andtables. The final decision is, therefore, a 
human one which requires time to develop, possesses inherent chances 
for error in the calculations, and, by its very nature, must be subject 
to the human bias and judgement of the decision-maker which may or 


may not adversely effect the decision. 


II, STATEMENT OF SPECIAL PURPOSE 


The advent of the computer has greatly facilitated scientific analysis 


and routine decision-making. To date the computer has not been applied 








to the search and rescue problem although the idea has been expressed 
and has been the subject of some feasibility studies. The application 
of computers in this field should greatly facilitate the analysis of the 
many factors involved and also reduce the time required to reach a 
decision. In view of the fact that numan lives are at stake in these 
operations, it is mandatory that the best possible decisions be made 
with the information and alternatives available, and these decisions 
must be made in the shortest possible time. The following problem has 


been selected with these thoughts in mind. 


III, THE PROBLEM 


It is the purpose of this study (1) to completely develop a computer 
flow diagram and associated computer program using the Control Data 
Corporation compiler Fortran 60 language (equivalent to Fortran II) for 
that portion of the "Distress" phase of search and rescue (SAR) operations 
relating to the determination of a datum point and associated search area 
Chapter Six; and (2) to investigate the possible development of those 
programs encountered during the research process which appear to offer 
improved techniques for analyzing and evaluating various portions of the 
problems encountered in the "Distress" phase. These potentially improved 
techniques will be submitted to the United States Coast Guard for 


evaluation. 








IV. DEFINITIONS OF KEY TERMS 


The major key term under consideration is the "Distress" phase. As 
to that phase of search and rescue operations wherein the estimated 
location of the distress incident and an associated surrounding area 
called the “search area” are determined and designated for any specified 
time. It also includes the assignment and dispatching of search units to 
the scene of the distress incident. However, for the purpose of this 
study no consideration will be given to the assignment and dispatching 
of search units. 

The following terms directly associated with the "Distress" phase 


will be used: 


Distressed unit. The distressed unit is that craft whose survivors 
require rescue or assistance. rs 

The object of the search may consist of a disabled vessel of any 
type, a downed aircraft, a life raft or life boat, or any form of survival 
equipment which will support a survivor or survivors on the water's 
surface, or, in some cases, the bodies of those who did not survive but 


which may be floating. 


Rescue Coordination Center. The Rescue Coordination Center (RCC) 


is that operational office responsible for the initiation, coordination and 








overall control of search and rescue operations for the specific areas 
under its jurisdiction. For Pacific operations this is controlled by 
Commander, Western Area, and for Atlantic operations it is controlled by 
Commander, Eastern Area. The specific individual responsible for the 
yarious decisions required on any given distress case is referred to as 


the RCC Controller. 


On-Scene Commander. This is the commander of a search unit at 
or near the scene of the search who has been designated by the RCC 
Controller as the commander responsible for controlling and coordinating 
the efforts of the various units participating in the search while they are 


Operating at the scene. 


search unit, Any unit which is actively participating in the search 
while under some form of control of the Rescue Coordination Center (RCC) 
or On-Scene Commander (OSC) is designated as a search unit. This may 


be an aircraft or surface vessel, commerical or military. 


Distress incident position. This is the geographical position of the 
distressed unit at the time that the distress ocurred. The accuracy with 
which this position is determined will depend upon the specific nature of 
the distress and the time and facilities available for determining the 


position. 








Datum. The datum point is that geographical position which represents 
the most likely location of the survivors at any specified time. Generally, 
the initial datum will correspond with the "distress incident position. " 

At any time after the distress has occurred, the datum point will have 
moved in response to the forces of Nature and the efforts of the survivors. 
The datum point is considered to be in the center of any designated search 
area. It serves as a reference point for determining the limits of a 


designated search area. 


Search area. That area in which the survivors are believed to be 
located and which can be covered with a reasonable degree of thoroughness 
by the search units is defined as the search area. It is bounded by limits 
established by the search radius which extends outward from the datum 


point. 


Search radius. The search radius, (R), is the distance (in nautical 
miles) from the datum point to the edge of the search area. Itisa 
function of the number of the search, for example=--the second search, 
and the total probable error of position. The number of the search determines 
the safety factor to be applied. The total probable error of position is a 
function of several other factors which are defined below. In general, if 
the search radius is R, the safety factor is SF, and total probable error 
of position is c, then 


R=SFxc 





Total probable arror of position, The total probable error of position, 


(c), is that distance (in nautical miles) which represents the best 
estimate of position error when the following factors are taken into 
consideration: 

a. Initial position error of the distressed unit, (X). 

b. Navigational error of the search unit(s), (Y). 

c. Life raft drift error, (de). 
Total probable error of position is a scalar quantity which is applied to 
vectors extending in all directions from the datum point. It is represented 
by the equation: 


c= Vdo2 +x? + Y2 


A review of the development of this equation is contained in Chapter Il. 


Initial position error, (X). This is an estimate of the distance (in 


nautical miles) by which the distress incident position may be in error. 
It is a function of the method of navigation used to fix the position of the 


distressed craft. 


Navigational error of the search unit, (Y). The navigation used in 


guiding the search units to the distress incident position or datum point 
and throughout the conduction of the search will be subject to varying 
degrees of error depending upon the method of navigation used, This 
error (in qauan miles) must be considered when determining the total 
probable error of position, c, and the resulting search radius. 


9 








Life raft drift error, Ge). Once survivors are free of the distressed unit 
and floating upon the surface, they will be subjected to both surface currents 
and surface winds. The interaction of these forces, the nature of the 
survival craft--that is, whether it is a life raft, life boat or life jacket-- 
and the time span over which they react, all contribute to a resultant 
drift vector, D. This drift vector may have maximum and minimum values 
depending upon the variations in the forces of Nature during the period of 
interaction and the availability of accurate information about the weather 
conditions at the scene of the distress incident. This variation introduces 
the concept of life raft drift error, (de), which will contribute to the total 


probable error of position. 


Drift. The movement of the survivors (datum) upon the surface is 
represented by the vector, D. This is the vector which is the sum of the 
average sea current vector, the local.wind current vector, and the leeway 
vector. Although the average sea current is assumed to have a constant 
value, the latter two vectors may vary thus contributing to maximum and 


minimum values for drift, D. 


Drift (max). Drift (max), Dy»sy, is that vector which represents the 
greatest displacement possible for a given datum during a specified 


period of time, At. 


Drift (min). Drift (min), D,j;,, is that vector which represents the 
minimum displacement possible for a given datum during a specified 


period of time, At. 
10 
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Average sea current. The average sea current is that vector described 
by the set and drift of the prevailing ocean currents for the area of the 
distress at the time of the incident. Itis determined by extracting the 
applicable values from the appropriate charts in the Atlas of Surface 
Currents* for the area to be searched. The average sea current vector is 


only considered if the Atlas of Surface Currents indicates that the reported 





current has a steadiness equal to or greater than 35% of the time. 


Local wind current. The direction and magnitude of the average sea 
current may be altered by the presence of local surface winds. The 
magnitude of the current generated by the local surface winds isa 
function of the force of the winds, their steadiness both in magnitude and 
direction, the size of the area over which they have been blowing (fetch), 


the length of time that they have been blowing, and the influence of the 


Ooriolis force. 


4Atlas of Surface Currents (Washington: U. S. Navy Hydrographic 
Office, 1953). This Atlas consists of a series of different atlases for 
various ocean areas of the world. For any given area the charts are 
furnished for each month of the year. These charts are divided into grids 
for each degree of latitude and longitude. An arrow and numbers placed 
in each one-degree quadrangle shows the mean direction of force of the 
surface current in that quadrangle for that particular month under average 
normal conditions and the number of observations reported for that 
quadrangle. A current rose shows the frequency distribution (steadiness) 
and average drift rate for the eight cardinal and intercardinal points for 
larger quadrangles. 
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Leeway. In addition to generating local wind currents, the local 
surface winds also act directly upon that portion of the survivors or their 
survival craft which is above the surface. The resulting displacement is 


called "Leeway. " 


Parachute drift. If the distress incident involves an aircraft distress 
which has required that the survivors bail out, consideration must be given 
to their free fall and parachute drift. The distress incident position will 
be the position of the aircraft at the time of bail-out, and the initial 
datum will be that position where the parachute would have landed after 
considering the altitude at which the aircraft was flying and the prevailing 


winds at different altitudes from the point of bail-out to the surface. 
IV. ASSUMPTIONS 


In approaching the first part of the problem it has been assumed that 
the procedures, graphs, tables and guidelines outlined in the National 
Search and Rescue Manual are based on the best information and theories 
available for the ep ctica! methods presently used by the United States 
Coast Guard in search and rescue operations. No further efforts have 
been directed toward verifying or cdisproving the existing doctrine. 
Instead, efforts have been concentrated on the development of - 
operational computer program which will facilitate the application of the 


existing doctrine. 


13 








In developing the computer program it was assumed that the various 
position and navigational errors could be approximated by normal 
distribution curves. Furthermore, where a vector produced two possible 
positions (see Figure 2, page 11), that is a maximum and a minimum 
position, it was assumed that the most likely position would be at a point 


midway between the two positions lying ona line connecting the two points. 


V. LIMITATIONS 


Certain self-imposed limitations were made in considering the 
problem, First, the problem was limited to the type where the time that 
the distress occurred was known, the nature of the distress and the type 
of units involved were known, and an estimated distress incident position 
was known. Furthermore, the estimated distress incident position was 
placed sufficiently far offshore to eliminate the effect of tidal currents, 
rotary currents resulting from coastal configuration, and currents caused 
by rivers flowing into the sea. 

Additional limitations precluded the development of the parachute 
drift program. Recent research into this area has resulted in the development 
of a computer program to handle this problem. Where an aircraft is 
involved as the distressed unit, the program is limited to the type of 
distress which involves ditching the aircraft into the ocean. 

No attempt was made to develop a program involving a small craft 


as the distressed unit since such a craft is unlikely to venture far enough 
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away from the coast to fall within the parameters adopted in limiting the 
distress incident position. 

Testing of the model was done by generation of hypothetical cases, 
manually computing a solution using the techniques and procedures 
these results with those produced on the computer. Historical data of 
actual distress cases were not available at the time that this project 
was completed. 

The program was developed for operations in the Atlantic Ocean, 
north of the equator only. With minor modifications, the program can be 
adapted for operations s@uth of the equator and in the Pacific Ocean, both 
east and west of the International Date Line. 

One final limitation on the study relates to the method used in 
measuring distances. All distances are treated as rhumb line distances 
based on Mercator Sailings rather than Great Circle distances. Since 
the magnitude of the distances involved is small, the error resulting in 


using Mercator distances will be neglected. 
VI. RESEARCH SIGNIFICANCE 


The development of an operational computer program to assist the 
decision-makers responsible for coordinating and conducting search and 
rescue operations will constitute a major step forward in improving Coast 


Guard operations. By eliminating the complex and tedious manual 
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calculations and plotting of information, it will free the Rescue 
Coordination Center Controller and his team to concentrate more fully on 
those phases of the operation which require the all important:tasks of 
human decision-making=-the final value judgments. 

The development of this program opens new horizons for the use of 
computers in expanding search and rescue operations by presenting a 
workable program which may be adapted to the existing Atlantic Automated 
Merchant Vessel Reporting System (AMVER) and the newly developing 
Pacific AMVER system. 

The most important significance of such a program is that it should 
lead to improved effectiveness in that mission of the United States Coast 
Guard which is directed to the preservation and safety of the lives and 


property of all those who travel over and upon the seas. 
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CHAPTER II 


REVIEW OF THE LITERATURE 


Studies directed toward the problem of search and rescue have been 
limited primarily to those conducted by the various Armed Forces and 
those specific agencies of other countries that have been responsible 
for the execution of this type of work. The United States Coast Guard 
has been the leader in this field both in the United States and in 
international efforts to develop the best techniques possible, Accordingly, 
the literature dealing with this problem has been produced by the combined 
efforts of the various services and agencies and certain individuals 
directly connected with these operations. 

Since the problem of search and rescue is amenable to study by the 
various phases involved, this review of the literature will be limited to 
that phase dealing with the determination of the best estimated position 
of the survivors at any given time having started with an initial estimated 
position, and then the determination of a reasonable area for searching 
around this estimated position, Furthermore, this review will be focused 


on two specific areas. The first will be devoted to the doctrine prescribed 


ype ee 


of Search Areas." It is for a portion of this Chapter that the computer 


program has been developed. 
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The second area deals with the literature which sheds new light on 
certain aspects of the problem which may tend to improve the existing 
doctrine as currently outlined in the National Search and Rescue Manual. 
Therefore, all literature reviewed will be that which has been developed 


reflecting the latest theories and ideas since the development of the 


National Search and Rescue Manual in 1961. 


SS ee 


I. DOCTRINE ON THE DETERMINATION OF THE SEARCH AREA 


ACCORDING TO THE NATIONAL SEARCH AND RESCUE MANUAL 


Before search and rescue units are dispatched to the scene of a 
distress and before the type of pattern of search to be used is determined, 
an estimated position of the survivors must be established, and an area 
sufficiently large to insure inclusion of the survivors must be selected. 
This search area is a function of the various parameters of the incident 
and the associated position errors. There are two possible situations 
confronting the decision-maker. First, the initial distress position is 
known with varying degrees of accuracy; and the second, the initial 
position of distress is completely unknown. Only the first case will be 
covered in this paper. 

There are many factors associated with determining the initial distress 
position and its related position errors. The following are some of the 


questions which must be answered by the search and rescue (SAR) 
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coordinator (RCC Controller). What type of craft is in distress? What 
is the nature of the distress? What type of navigational information was 
used to determine the initial distress position? What type of survival 
craft is likely to be employed by the survivors? What are the existing 
weather conditions at the scene of distress? What has the weather been 
like at the scene for the preceding twenty-four hours? How long will it 
be before a SAR unit can reach the scene? Where is the most likely 
position (datum) of the survivors when the SAR units reach the scene? 
What are the various drift factors which must be considered in this 
particular case? What is the expected time of survival forthe survivors 
considering sea water temperature and surface weather conditions at the 
scene? 

first divides the factors affecting the position into three broad classes. 
The first of these is called the "Initial Error in Position"» and is a function 
of both the type of craft involved and the method of navigation used. If 
the reported initial position of the distress incident is based upon an 


accurate navigational fix, the following position errors are assumed:® 


5 
National Search and Rescue Manual, op. cit., p. 6-2. 


Slbid., p. 6-7. 
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Radius of 


Type of distress craft Assumed Position Error 
Ship Fix Five nautical miles 
Aircraft Fix Ten nautical miles 
Small Craft Fix Fifteen nautical miles 


If the initial position of the distress incident was determined by monitor 
stations on radar nets or on military or commerical radio direction-finding 
nets, these various nets will classify their fix. These classifications 
have standard specified position errors. 

If a fix could not be obtained at the time of the incident, or if the 
reported position was not specifically reported as a fix, it is assumed that 
the position is based on dead reckoning (DR). In this case the following 
pasition errors are assumed:/ 

Type of distress craft Radius of Assumed Position Error 

Five percent of the distance 
Ship DR traveled since the last fix, 

plus the error of that fix. 

Ten percent of the distance 
Aircraft DR traveled since the last fix, 

plus the error of that fix. 

Fifteen percent of the distance 
Small Craft DR traveled since the last fix, 

plus the error of that fix. 

The second factor affecting the distress position is that of parachute 


drift where the particular distress incident involves an aircraft from which 


the survivors have bailed out. The position reported as the initial 


Tid, 
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position of the distress incident will not be the same position as that of 
the survivors when they land in the water. A table is provided to determine 
how far downwind the landing will be from the initial reported position of 
distress or bail-out paint. The entering arguments for this table are 
altitude of parachute opening and the wind velocity in knots which 
represents the average wind between the winds at the altitude at which the 
parachute opens and the surface. 8 

The third factor, and the one most subject to question, is that of 
"Survival Craft Drift."9 Once the survivors are adrift, they are subjected 
to the various elements of Nature. In the open sea the simplifying 
assumption has been made that the resultant movement of the survivors 
is the consequence of three forces acting upon them. These three forces 
influence the survivors motion or travel independently. They are (1) the 
force due to the average sea current which is a function of the major 
ocean current circulation pattern; (2) the force due to local wind currents 
which have been generated by the local winds in the area for the twenty- 
four hours preceding the distress incident; and (3) the leeway of the 
survival craft as it responds to the local winds at and following the time 


of the incident. 


8Ibid., p. 6-3. 
ding. 
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The average sea current is determined by consulting the appropriate 


Atlas of Surface Currents for the area and month involved. If the average 








sea current is indicated as having a steadiness (frequency of reports) of 
35% or greater, it is considered as a relevant force. Otherwise, it is 
disregarded. 

The local wind current is treated as follows: 

a. Obtain an estimate of the local wind vector for the preceding 
twenty-four hour period. If the vector has not varied, assume that the 
local wind current has diverged from the wind vector by thirty degrees to 
the right for latitudes north of 10°N, thirty degrees to the left for latitudes 
south of 10°S, or is directly downwind for latitudes 10°N to 109S. This 
divergence reflects the effect of the Coriolis force. 

b. If the wind vector has varied during the twenty-fottrhour 
period, determine the mean wind for four hour periods during the twenty- 
four hours. Then use the mean wind to determine the direction.of the 
local wind current. If variations have been more than 45° or ten Knots, 
intervals smaller than four hours should be used, with the final local wind 
current direction becoming the sum of the vectors:forall the time intervals. 

c. The magnitude of the local wind current is taken froma table 
with an entering argument of the local wind speed in knots. 

d. The assumption is made that for the first six to twelve hours that 


a wind is blowing, the local wind current generated will be downwind. The 
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divergence indicated ina. above applies only after the wind has been 
blowing for more than twelve hours. 

The component of leeway is equally complex. For survivors supported 
by life jackets or floating free on the surface, the sail area exposed to the 
winds is considered to be negligible as compared to the submerged surface 
area which will be subjected to the surface currents. For those survivors 
supported by life rafts, the leeway will be directly downwind. The 
magnitude of this leeway will be dependent upon whether or not the 
survivors have streamed a drogue. If the survivors are ina life boat, the 
magnitude will also be dependent upon the absence or presence of a drogue. 
In addition, the drift of a life boat may vary as much as forty degrees on 
either side of a downwind direction depending upon the configuration and 
loading of the life boat. 

If the winds have been fluctuating, the leeway is determined for four- 
hour periods and then added vectorially to determine the final component. 
the magnitude of a typical life raft leeway (the type is not indicated) for 
conditions with and without a drogue. 10 

Whenever uncertainty exists as to the length of time adrift, the 
magnitude of the leeway, or the accuracy of the reported conditions of 


wind and current, leeway is approached using minimax plotting techniques. 1] 


10Ibid., p. 6-4, 
» ; : 1) Se 
Mibid., p. 6-5. 
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Under these conditions the RCC Controller must determine the minimum 
and the maximum probable values of wind and current. From this information 
a plot is constructed to reflect the total minimum possible drift and the 
total maximum possible drift. The "most probable location of the target 
is the position midway between the minimum and maximum drift positions. " l2 

If, however, the wind conditions have been fairly steady, an average 
value may be used. This will lead to a final value of the drift vector which 
is then considered to be an average drift vector. Both of these concepts 
were illustrated in Figures 1 and 2,Chapter I, page ll. 

The status of the problem at this point finds the RCC Controller with 
an estimated position of the incident based on information received in the 
"Alert" phase. He also can determine the most probable location of the 
Survivors after applying the drift factors. The next step is to evaluate 
the total probable error associated ne the latter position. 

The total probable error of position is obtained by considering the 
initial position error, the navigational error of the search unit(s), and the 
drift error associated with the determination of the drift of the survivors. 

The first of these, initial position error, is determined in the manner 
described earlier in this Chapter. The navigational error of the search 
unit(s) is handled in a similar manner using an identical table for 


determining the value of the error. The values for these two errors are 
12 ibid. 
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represented by "X" and "Y" respectively. 

The drift error, de, can be determined one of three ways. If the 
total drift was determined as an average drift, the drift error is estimated 
to be equal to the total drift divided by eight. 13 If minimum and maximum 
drifts have been computed, a table is provided for determining the drift 
error using the minimum and maximum drift values as entering arguments. 

The third method of determining drift error is used when the maximum 
and minimum drifts have been plotted. Using the values for the maximum 
and minimum drift, the respective drift errors are determined by dividing 
each value by eight. These values for de are then plotted as the radii of 
Circles centered on their respective drift positions. The final drift error, 
referred to as the minimax drift error, is found graphically to be the 
radius of a circle centered on the line connecting the maximum and 
minimum drift positions and which is tangent to the outer limits of the 
two circles generated by the maximum and minimum drift errors. See 
Figure 3. 

The total probable error of position can now be determined. This 


error, Cc, is determined by the equation: 


c= Vd? + X44 y2 
l3ipid., p. 6-7. 
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Although the National Search and Rescue Manual does not present the 


development of this equation from the theoretical background, a separate 


14 


study~* conducted in 1964 derived it mathematically. The equation was 


derived using the Law of Cosines and the root-mean-square of maxima 


and minima treating the three error vectors, two at a time. 
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FIGURE 3 


GRAPHICAL SOLUTION OF DRIFT ERROR (MINIMAX) ~ 


Furthermore, if the values for X, Y and d, are considered to represent 
the standard deviation for a normal distribution associated with each of 


these errors, the variance would be X2, Y2 and do respectively. Then c 


ye obert L. Armacost and Norman T..Saunders, "Computer Program 


of USCG Search and Rescue Procedures" (unpublished Bachelor's research 
paper, The United States Coast Guard Academy, New London, 1964). 
unnumbered. 
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would represent the standard deviation of a normal distribution representing 
the sum of these errors. The 1964 study of Armacost and Saunders also 
reached this conclusion. !° 

The American Practical Navigator!6 devotes a chapter to the 
consideration and study of navigational errors. Among the errors included 
in this discussion are constant errors, random errors, periodic errors, and 
systematic errors. 

Then attention is directed to the combination of these errors in any 
given line of position or navigational fix. The distribution of these various 
errors and the application of probability of error is applied. The resulting 
conclusion is, "If a number of random errors are combined, the result 
tends to follow a normal curve regardless of the shape of the individual 
errors....."1/ Furthermore, applying statistical theory the authors state, 
_",..the standard deviation of all the errors is found by squaring each 
individyal error, adding the results, and taking the square root of the sum." 18 
They further point out that a given line of position, as finally determined 


and plotted, "...may include 30 errors or more. " 19 


1SIpid, 

16Nathaniel Bowditch and others, American Practical Navigator, 
H.O. Pub. No. 9, (1958 edition, Washington: U.S. Navy Hydrographic 
Office, 1958), pp. 678-688. 

‘7ibid,, p. 683. 

18Ipid, 

1Sipid., pe 684. 
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Although Dutton, in his book, Navigation and Nautical Astronomy, 
Goes not comment directly upon the distribution of navigational errors, 
he does discuss celestial fixes for surface craft and the magnitude of the 
errors encountered. He states that, "...under average conditions an 
experienced navigator can obtain a fix that is seldom in error by more 
thana mile."29 For celestial navigation in aircraft, "Under good 
conditions...a fix accurate to 5 to 10 miles is considered good..." but 
"...individual sights taken by experienced navigators may be more than 
100 miles in error in extreme cases, "21 

In considering loran fixes he states: 

"In general, ground waves produce a line of position 

accurate to 1.5 miles or better over 80% of the area 

covered by this stations. Sky waves produce a maximum 

error of 5 to 7 miles for 80% of the coverage area. The 

accuracy of a position is increased, of course, if three 

or more lines of position are used and the various lines 

weighted by their position relative to the stations, type 

of wave, angle of cutting other lines, etc. nd 

Having determined the total probable error of position, the next step 
is to use this information in determining the limits of the area to be 
searched. Initially, the search area must be sufficiently large to insure 
that the total probable error of position is included. The theoretically 
ideal search area would be a circle centered at datum and with a search 


20Renjamin Dutton, Navigation and Nautical Astronomy (tenth 
edition, Annapolis: United States Naval Institute, 1951), p. 332. 


éd tid. , Pp. 337, 
22ipid., pe 228. 
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radius equal to or greater than the total probable error of position, c. 

The search area is expanded for each subsequent search following 
an unsuccessful search by increasing the search radius, R. The number 
of searches conducted is not limited to a specific number since adverse 
weather conditions at the scene may reduce the effectiveness of each 
sweep. Such conditions would, therefore, necessitate additional 
searches over and above the number used under good search conditions. 

Each succeeding search will be centered on the datum point at the 
time the particular search is to be started. This means that the area where 
the survivors are most probably locatgfd, namely datum, will be searched 
more frequently than the outer limits of the search area. The outer limits 
are extended for each succeeding search. The search radius, which 
determines the outer limits, is found by multiplying the total probable 
error of position, c, by a safety factor. The safety factor is increased 
for each search to reflect the increased total probable error of position 
which has resulted from increased drift and the associated increased drift 


error. The values used in determining the search radius are as follows:23 


search Safety Factor Search Radius 
] el lL loxge 
Z LG 16xCc 
3 AAG, 2.0xc 
4 Ze 2.3 XC 
3 2.9 PARES >< (6: 


23National Search and Rescue Manual, op. cit., p. 6-10. 
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The time required to complete a search will be a function of many 
things, such as the number of search units participating, the speeds 
used by the search units, the type of pattern used in the search, and 
the size of the area to be searched. If several hours elapse in conducting 
a search, it may be necessary to completely recompute the total probable 
error of position, c, ,to reflect the increased drift error. 

As indicated earlier, the theoretical search area would be a circle. 
However, to attempt to search while navigating within the limits ofa 
circle is not particularly suitable for many of the different types of 
search patterns commonly used. Therefore, the search area is usually 
squared off to provide a rectrangular search area with a side of length 
2R and a total area equal to 4R2, This search area provides a more 
suitably defined area for navigational purposes. 

Figures 4 and 5 show the concept of the "repgated expansion" search 
concept for a stationary datum and a moving datum. This concept calls 
for five searches following consecutively using the search radii indicated 
previously. This provides for the greatest concentration of search efforts 
in the area where the survivors are most probably located and minimum 


efforts in the areas of least probable location. 24 


24ibid., p. 6-12 
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FIGURE 4 


REPRHATED EXPANSION SEARCH AREAS FOR STATIONARY DATUM POINT 


5th search 
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REPEATED EXPANSION SEARCH AREAS FOR MOVING DATUM POINT 





Having determined the datum point and search area, the ROC Controller 
must then determine what search units are available, how long it will take 
them to reach the scene, what type of search pattern is most appropriate, 
and how much time is available for searching (limiting factors include 
aircraft fuel and hours of daylight remaining before dark). Considering 
these factors, the ROC Controller must make a value judgment of the 
search area calculations and perhaps modify the area by altering the size 
of the search radius or the track spacing of the search units or both. This 
is known as a "feasibility check."2% This is a decision which must be 
made on the basis of the specific situation whicn exists for each individual 


distress incident. 
II. STUDIES DEALING WITH PARACHUTE DRIFT AND OCEAN CURRENT DRIFT 
Parachute drift. 


In 1964, while at the United States Coast Guard Academy, Cadets, 
First Class, Armacost and Saunders, conducted a study29 of the general 
problems associated with search and rescue operations. A major portion 
of their efforts was directed in examining the problem of free-fall and 
parachute drift of survivors who were required to bail out of a distressed 


aircraft, 


2°Ibid., p. 7-11. 
20r raacost and Saunders, Op. Cit. 
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Their examination of the table provided in the National Search and 
Rescue Manual lead them to believe that the values contained therein were 
related linearly. They also concluded that the table was based on a model 
where the survivor's initial velocity at time of bail-out was zero in both 
the horizontal and vertical planes. Also, the table utilizes an average 
wind between bail-out altitude and the surface. It makes no provision for 
the wide variance in both wind velocity and direction which may be 
encountered at different altitudes. 

Accordingly they set out to derive equations which would relate 
the force of gravity, effect of winds in a horizontal direction, motion due 
to the velocity of the aircraft, and air resistance. They did this by 
developing equations for the vertical component of free fall, the horizontal 
component of motion due to wind, the horizontal component of motion due 
to an initial velocity equal to that of the speed of the aircraft, the effect 
of air resistance or drag as a function of air density which in turnisa 
function of barometric pressure and temperature, the relationship of 
kinematic viscosity of air, air density and altitude, and finally, the 
evaluation of the various constants appearing in the various equations. 
Having completed this extensive development, they then wrote a basic 
flow diagram for their program. That portion of their work dealing with 
the parachute drift problem has been reproduced in its entirety in 


Appendix C. 
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In 1965, Cadet, First Class, Darvy M. Cohan, also of the United 
States Coast Guard Academy, continued the work started by Armacost and 
Saunders. Cohan refined and developed a computer program2/ incorporating 
the results of Armacost and Saunders; but, in addition, he investigated the 
effects of changing wind direction and velocity using three different 
altitudes: (1) surface winds, (2) medium altitude winds, and (3) high 
altitude winds. These winds were selected since information on their 
magnitudes is usually available from various agencies providing 
meteorological data. 

In addition, he considered the effects of three possible parachute 
opening conditions which included the two possible extremes, namely (1) 
that the parachute opened immediately upon leaving the aircraft, or (2) 
the parachute did not open at all, and the possibility of opening 
automatically at 14,000 feet. This latter altitude is the one at which most 
military parachutes are designed to open using a barometric automatic 
release mechanism. 

Cohan employed the Euler method of integration in his program. 
He tested his program on an IBM 1620 computer and found that the running 
time was slightly less than ten minutes. Accordingly, he recommended 

27 Dar vy M. Cohan, "Free Fall and Parachute Drift" (unpublished 
Bachelor's research paper, The United States Coast Guard Academy, New 
London, 1965). 
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that future work in this area should include consideration of utilizing more 
sophisticated techniques such as the Runge-Kutta method. Furthermore, he 
pointed up the weakness that exists in that no consideration has been given 
to the effect of turbulent weather conditions with the associated vertical 
shear, updrafts, and downdrafts which may have contributed to precipitation 
of the distress itseli. 

Appendix D contains the basic equations, flow diagrams, and computer 


program for the IBM 1620 which Cohan developed. 


Ocean current drift. 


As military and commercial maritime operations continue to grow, 
the need for more and more accurate information dealing with surface (and 
sub-surface) currents~-their prediction and behavior~--increases in 
importance. Oceanographers and meteorologists have long been studying 
this problem. Many theories have been advanced; however, obtaining 
sufficient experimental data to test these theories has proved to be a 
major obstacle. 

The Fleet Numerical Weather Facility (FNWF) at Monterey, 
Califomia, has recently developed a computerized system for producing 
synoptic analyses of surface currents based on the distribution of sea 


surface temperatures and the climatological temperature field at a depth 
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of 200 meters. Commander W. E. Hubert, USN, presented a paper28 at 
the Institute of Navigation in San Francisco in December 1964 is which 
he discussed the general method used in the Fleet Numerical Weather 
Facility’s program. This program is designed to evaluate the transport 
due to the permanent flow component and the component due to wind and 
waves thus producing the overall pattern of surface currents. Commander 
Hubert’s paper appears in Appendix E, 

Taivo Laevastu has conducted an extensive review of the materials 
dealing with surface currents in connection with his work at the Fleet 
Numerical Weather Facility and his preparation of the book Oceanographic 
Forecasting for the Navy Oceanographic Office. In Chapter Ten, Doctor 
Laevastu examines in detail the various aspects of "Prediction of Surface 

_ Currents. 29 Among the factors investigated, he discusses the evaluation 
of permanent flow, wind currents, mass transport by waves, inertia currents, 


tidal and "hydraulic currents, " currents caused by the change in atmospheric 


pressure, and other current components. 30 He further discusses the 


28 E. Hubert, “Computer Produced Synoptic Analyses of Surface 
Currents and Their Application for Navigation" (paper read at the 1964 
National Marine Navigation Meeting, Institute of Navigation, San 
Francisco, December 7-8, 1964). 


2373 ivo Laevastu, "Prediction of Surface Currents, " Oceanographic 
Forecasting (a book presently under preparation for U.S. Naval 
Oceanographic Office under Contract Number N62306-1129). 

30ipid, , Ppp. 10-1 = 10=26. 
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evaluation of change of current speed and direction with depth, the 
influence of water depth, configuration of coast and islands, nearshore 
_, circulation, current boundaries and eddies, and deep currents. 31 
Further studies relating the work of the Fleet Numerical Weather 
Facility and the problems associated with surface weather phenomena 


are found in their Technical Notes Numbers Fiyes2 and Bight?3 published 








in February 1965. These Technical Notes discuss the sea-air interactions 





on.a synoptic scale and the analyses and forecasting of sea surface 
_ temperatures, both of which are directly related to the paper presented 
by Commander Hubert in December 1964 and referred to earlier. 

The U. S. Naval Oceanographic Office has been vitally interested 
in the problems associated with ocean current prediction. Although 
military operations involving anti-submarine warfare have given impetus 
to. studies in this area, the problems associated with prediction of the 

_ location of freely drifting mines and rafts and the associated areas for 
searching for these objects have also been under investigation. 


Donald A. Burns of the U. S. Naval Oceanographic Office has 


ee. oak 10-24 2 onsey 


Say, E. Hubert, "U. S. Fleet Numerical Weather Facility Activities 
Relating to Sea~Air Interactions on a Synoptic Scale” (Technical Note 
Number Five. Monterey: Fleet Numerical Wé@ather Facility, 1965). 


33p_ M. Wolff, L. P. Carstensen and T. Laevastu, "Analyses and 


_ Forecasting of Sea Surface Temperature (SST)" (Technical Note Number 
Eight. Monterey: Fleet Numerical Weather Facility, 1965). 
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developed a Fortran Computer Program which defines a drift ellipse pattern 
for the free-floating object. 34 The essentials of this drift ellipse program 
are described inA Statistical Rose Program developed by Walter E. Yergen 
and published by the Oceanographic Office. 39 

Basically, Burns' program analyzes input data consisting of 
“grouped frequency current observations. "36 This analysis produces 
various parameters which are then plotted to outline that area in which 
recovery is most probable. This is accomplished by determining vector 
means, vector variances and covariances, determining the inclination of 
the principle axis of the distribution, and the standard deviations along 
these principle axis. Next, within specified limits for speed and direction 
classes, the frequencies for the normal surface are integrated. 37 This 


Drift Ellipse Program is found in Appendix F. 


III NOTE ON ADDITIONAL LITERATURE RELATING TO THE 
OVERALL SEARCH AND RESCUE PROBLEM 
In undertaking this project it was necessary that a thorough general 
34Donald A, Burns, "Instructions for Surface Drift Ellipse Program" 
(Washington: U. S. Naval Oceanographic Office, June 1962). 
3°Walter E. Yergen, A Statistical Rose Program, SP-64 (Washington: 


Evaluation Branch, Oceanographic Analysis Division, Marine Sciences 
Department, U. S. Naval Oceanographic Office, October 1962). 





3©purns, op. Cit., passim. 
37Yergen, op. cit., passim. 


38 








understanding of all aspects of the search and rescue problem be acquired 
before an attempt was made to sub-optimize that portion which was 
selected. Accordingly, the Bibliography at the end of this paper is 
annotated to reflect all the pertinent literature which was reviewed in 


the process of this study. 
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CHAPTER [I 


THE STUDY 


I. METHOD 


The method of research employed in this study was basically that of 
the development of a conceptual model designed for a computer solution 
expressing the doctrine outlined in the National Search and Rescue Manual 
for the determination of a search area. Through the application of 
mathematical techniques relationships of the numerous variables were 
developed. These relationships were then outlined in a series of computer 
-flow diagrams, for various subroutines, and the computer programs were 
written in Fortran 60 computer language. After each subroutine was 


thoroughly checked and tested, the final program was assembled, tested 


and evaluated using hypothetical distress cases. 


Il. PROCEDURE 


The problem of determining a Search area for any given distress 
incident involves three basic steps: 

1. Determine the datum point at any specified time. 

2. Determine the probable position errors associated with the datum 
point. 

3. Determine the limits of the area which will be included in the 
search, 
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Each of these steps was examined individually and in detail. The 
examinations were accomplished incrementally by determining what 
output was required, what were the inputs available or necessary to 
develop the desired output, and how were the inputs and output related. 


The limitations outlined in Chapter I were carefully observed. 
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THE BASIC PROCESS OF DETERMINING THE SEARCH AREA 


The following sections;of this chapter will treat each of the above areas 
in order. Program flow diagrams will reflect the general logical sequence 


which was employed although the terminology used will be broad enough 
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to permit the development of similar programs using computer languages 
other than Fortran 60. Standard symbols will be used throughout the 
diagrams. An index to these symbols will be found in Appendix A. The 
detailed Compliter program in Fortran 60, complete with definitions of the 


applicable variable names used, will be found in Appendix B. 


Determination of datum at any specified time, tj. 

When a distress case has entered into the "Distress" phase, the 
RCC Controller's first task is to determine the datum point. Since this 
datum point will be continuously moving as time passes, he must select 
the specified time, tj, for which datum is to be determined. The initial 
reported position of the distress incident will be datum tg where ty is 
the time at which the distress incident occurred. The time selected, tj, 
may be the estimated time of arrival of search units at the scene, or it 
may be the predicted position at any other time during the period of the 
case duration. It will be redetermined during each search and used as 
the center of subsequent searches. What then is the information needed 
by the ROC Controller to determine datum tj? 

The basic information needed will be the initial Aienoee incident 
position, datum to, or the previous datum if this is a redetermination 
of datum, the time of the previous datum, the time for which the new 


datum is being determined, and the drift vectors which have been acting 
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upon the previous datum during the time interval which will have elapsed. 


See Figure 7-A. 
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THE FACTORS ASSOCIATED WITH THE FORECASTED DATUM 


? 


The following specific information is necessary: 

1. The type of unit in distress, i.e., an aircraft, a ship, ora small 
craft (this latter portion was not developed in this study). 

2. The initial distress incident position. This position will be expressed 


in terms of latitude and longitude. 
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3. The time associated with the initial distress incident position in 
terms of the date, hour and minute expressed in Greenwich Mean Time 
(GMT). 

4. Based on the type of unit in distress and a knowledge of the type 
of survival equipment carried on that type of distressed unit, the ROC 
Controller must specify what type of survival craft is believed to be the 
object of the search. This information may have been furnished in the 
message reporting the distress incident. 

5. Weather information at the scene of the distress incident including: 

a. The present surface wind velocity expressed in knots and the 
. direction from which it is blowing in degrees from true North. 

b. The surface winds for the twenty-four hour period preceding 
the time for which the datum is being determined. The velocity and 
directions are needed and must be expressed in knots and degrees from 
true North. 

6. The average sea current set and drift expressed in knots and 


degrees from true North. This information is extracted from the Atlas of 





Surface Currents for the appropriate area and month. At the same time 
the steadiness of this current is also determined and noted. 

7. The specific time, tj, for which datum is to be determined. This 
will be expressed in terms of date, hour, and minute using GMT. 

Using these specific inputs, the drift components of leeway, local 


wind current, and average sea current are determined as follows: 
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Leeway. Since leeway is the motion of the survival craft resulting 
from the action of the surface wind upon the craft, the first factor to be 
considered is the type of survival craft believed to be employed. If the 
survivors are believed to be floating in the water, either unaided or with 
life jackets, the effect of the local surface winds is assumed negligible. 
In this case the leeway component becomes zero and the computations 
continue to the next part of the program. 

If, however, the survival craft is believed to be either a raft ora 
boat, a leeway component will exist and must be determined. Considering 
the case of a raft first, its motion will be directly downwind. Therefore, 
the input of the local surface wind direction plus 180 degrees will yield 
the direction of the leeway component. The rate at which the raft will 
drift is a function of the local surface wind velocity and also depends 
upon whether or not the survivors in the raft stream a drogue to retard 
their motion in leaving the site of the distress incident. Therefore, two 
possible drift rates must be considered; (1) the rate of drift with a drogue, 
and (2) the rate of drift without a drogue. This results in the determination of 
a maximum value and a minimum value of the leeway component. When the 
leeway components (in knots) are multiplied by the elapsed time, At, the 
total maximum and minimum leeway components are determined. 

The treatment in the case of a boat is similar to that of a raft. 


However, the direction of drift of a boat has been found to vary as much 
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as 40% on either side of the downwind direction. Therefore, the maximum 
and minimum leeway vectors must also include maximum and minimum 
vectors 409 on either side of the downwind direction to reflect this 
divergence. This provides a total of six leeway vectors of which the 

two downwind vectors are used in computing datum while all six vectors 
are considered in the determination of the drift error. 

In this specific study, the magnitude of the minimum leeway component 
was determined by one of two equations based on leeway curves given in 
the National Search and Rescue Manual. 38 Analysis of the curves showed 
that leeway for a raft with a drogue could be closely approximated by two 
Straight lines. The first line approximated leeway for wind velocities 


equal to or less than five knots: 


Leeway in n. miles = (wind velocity in knots x 0. 8) 
ice 24 


The second line approximated leeway for surface wind velocities greater 


than five Knots: 


Leeway in n. miles = (wind velocity in knots + 2.5) 
Ae: 1.87 x 24 


The magnitude of the maximum leeway component (no drogue streamed) 
was also determined from the curves. However, the leeway curve values 
had been stored in memory and were furnished directly as a function of the 
wind velocity. Figue 7-B is the flow chart for the program used in 


determining the leeway components. 


38 
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FLOW DIAGRAM FOR DETERMINING LEEWAY DRIFT COMPONENTS 
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Applying the elapsed time to the leeway magnitudes and then 


combining this product with leeway direction yields the leeway vectors. 


Local wind current. The local wind current is treated in the manner 





described in Chapter II, However, the mean wind direction and velocity 
are determined manually with the number of calculations involved depending 
upon the extent of the surface wind variations during the preceding twenty- 
four hour period. The mean wind for the preceding twenty-four period is 
then used as an input for the computer. 

The magnitude of the current generated by the local wind is a function 
of the mean wind velocity. Examination of the table relating the velocity 
of the wind current to the mean wind speed in the National Search and 
Rescue Manual suggested a linear relationship. 39 The data from the table 
were plotted and the following equation provided the relationship: 


Velocity of the local wind current = (mean wind/1. 18) 
(in knots) 24 


The latitude of the initial distress position (or previous datum) was 
compared to the 10°N latitude criterion to determine the a ereanee of the 
local wind current resulting from the Goriolis force. If the distress ingldent 
latitude was less than or equal to 10°N, the direction of the local wind 
current became the direction of the mean local wind + 180°. If the distress 


incident latitude was greater than 10°N, the direction of the local wind 


Stata. 
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current diverged 30° to the rignt of the downwind direction of the mean 
local wind. 

By applying the elapsed time to the local wind current velocity and 
then combining this with the direction of the local wind current, the 


final local wind current drift vector was determined. See Figure 7-C. 


Average sea current. This component is simply the values determined 


from the Atlas of Surface Currents for set and drift with the latter portion 





multiplied by the elapsed time. If the degree of steadiness was less than 
35%, this component is assigned a zero value by the RCC Controller. 


See Figure 7-D. 


Resultant total drift, The final resultant total drift vector for the 
elasped time period is the vector sum of the leeway components, the 
local wind current component, and the average sea current component. 
Since this involves maximum and minimum components, maximum and 
minimum total drift vectors will be generated. In accordance with the 
assumptions stated in Chapter I, the new datum will be found at a point 
midway between the maximum drift position and the minimum drift position 
found by using the downwind leeway vectors. See Figure 2 in Chapter I, 
page 11, for the general drift picture. Figure 8 shows the drift picture 
for a life raft, and Figure 9 shows the drift picture for a boat with the 


associated divergence of the leeway vectors. In the latter case the most 
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probable position of the survivors is assumed to lie midway between the 
maximum and minimum positions on the downwind vector. However, the 
extreme drift positions associated with the divergence of the leeway 
vectors serve as input data for a portion of the problem phase relating 
to position errors. 

In the specific program developed by this study the problem of 
vector addition was treated in two steps. The first step was to canvert the 
directions of the vectors, which were expressed in degrees measured 
clockwise from True North, to the standard Cartesian coordinate system 
used by the computer. Computer angles are measured in radians. The 
conversion was accomplished by the relationship: 

Bp = 27 - OR 
Where BR is the angle expressed in radians as used in the computer 
calculations, and Op is the direction of the component converted from 
degrees True into radians. 

The second step in treating the vector addition problem was to convert 
the various components of leeway, local wind current, and average sea 
current into their respective x and y components. The final resultant 
vector is obtained by summing the x components, summing the y components, 
and then applying the Pythagorean Theorem to find the resultant vector. 

Actually, the various drift vectors were broken into the x and’y 
components and carried as subscripted variables throughout the program 


and applied where needed. 
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Determination of total probable error of position. 

The total probable error of position is a function of three basic errors: 
(1) the initial error associated with the distress incident position, (2) the 
navigation error associated with the search units, and (3) the error 


associated with the determination of the drift components. See Figure 
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COMPONENTS OF TOTAL PROBABLE ERROR OF POSITION 


The following specific information is needed to determine the values 


of these errors: 
1. The type of unit in distress. 
2. The type of search unit which will be dispatched... 
3. The classification of the method of navigation employed to determine 


the initial distress incident position. If the method of navigation is 
unclassified or is classified as dead reckoning (DR), the last known fix 


of the distress unit must be ascertained. 


o4 





4. The classification of the method of navigation which will be 
employed by the search unit in proceeding to the scene of the distress 
incident. 

5. The present position of the search unit must be determined. This 
present position of the search unit is treated as a fix. It must be expressed 
in terms of latitude and longitude. 

6. The magnitude of the maximum drift vector. 

7. The magnitude of the minimum drift vector. 

Applying these specific inputs, the three types of errors are determined 


as follows: 


Initial position error of the distress incident. This error is strictly 
a function of the type of craftin distress and the method of navigation used 
to determine the initial position of the distress incident. In a general 
program the procedure is to test for the type of craft in distress which then 
determines which portion of the program is to be used. The next step is 
to determine the classification of the method used in determining the 
reported incident position. Using the values for initial position error 
cited in Chapter II for the appropriate navigational method produces the 
desired initial position error, X. 

If the classification of the navigational method is that of a position 


based upon radio direction-finder bearings (DF), the type of DF system 


used must be further identified. The two types of DF systems commonly 
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used are the commercial-military DF service nets and the Federal 
Communications Commission (FCC) nets. Each of these nets will furnish 

a classification of their position fix which indicates the error associated with 
tne position. These classifications, which are used in the program, are 


given below. 49 


Fix Classification Associated Position Error As A Radius In Nautical Miles 


EEO «© ee ee ee a ee 





DF Service Net FCC-DF Net 
ClassA Within 5 20 or less 
Class B Within 20 40 or less 
Class C Within 50 60 or less 
Class D -~--- more than 60 


Since the DF classification system is independent of the type of 
craft involved, one portion of the computer program handling this type of 
position classification serves for all other branches of the system 
regardless of the type of craft invoived. 

Originally the specific program developed by this study followed the 
general program shown in Figure 10-B. The various navigational 
classifications were assigned different code numbers or combinations of 


codes. These codes were then punched into the input data cards. Since 


40tbid., p. 5-7. 
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it was obvious that the operator in the ROC would have to put in these 
classification code numbers in one form or another, it was decided to 
change the code numbers to the actual values of the error for the methods 
which were classified as a Fix or DF Fix. This limited the program 
calculations to those involved with DR classifications. This will probably 
be the classification most generally encountered in actual practice. 

In the case of a DR classification, the error is equal to the error 
associated with the last known fix plus a certain percentage of the distance 
traveled from that last fix to the reported DR position. The percentage of 
the distance used is dependent upon the type of craft involved. The 
problem encountered in this determination is ascertaining the distance 
traveled between those two points. 

There are various methods available for determining distances between 
points on the Earth's surface. The method selected for this program was 
that of Mercator Sailing. The only inputs needed for the computer were 
the latitudes and longitudes of the two points. The basic equations used 
in solving this problem were those developed by Bowditch in his book, 
American Practical Navigator. 4] They are: 


i 


- we Dis 
De Gosc ! and tanOG = “mn 


a 


a ccaabinela, Op.Cltis, perzcr. 
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where: D = the distance in nautical miles between 1 and 2. 


yk 


DLo = the difference in longitudes, 1 and Ao: in minutes. 


Il 


the difference in latitudes, L) and L,, in minutes. 


mM 


the difference in meridional parts, M, and M,, (where 
meridional parts are a measure of the arc length ofa 
meridian between a given latitude, Lj, and the equator 
measured on a Mercator chart, expressed in units of 

1' of longitude at the equator). 4 


M, = the number of meridional parts between the equator and 
ayinsuvetes Uae 


The values of Mj were obtained by using the first two terms of the 
equation: 43 
Mj = 7915.704468 log tan(45°9 + Lj/2) - 23. 268932sinL; 


- 0.052500sin3L; - 0.000213sin°L; - 


of the search unit. This error follows the same 


Navigational error 





identical pattern of solution as that utilized in the determination of the, 
initial position error, X. The general flow diagram is shown in Figure 
10-C,. Only the position of the last known fix of the unit concerned will 


vary. 


42tbid., p. 936. 


43Ibid., p. 1187. "The constants used...are based upon Clarke's 
spheroid of 1866...the standard reference spheroid used for charting 
North America." In the computer program only the first six significant 
digits were used for each constant. 
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Drift error. The methods available for the determination of drift 
error, d,, were outlined in Chapter II. The computer solution for de 
was achieved by using the vector summation of the various components 
of drift as discussed earlier. These vector sums produced maximum and 
minimum drift vectors. In the case of a raft there are two vectors, one 
maximum and one minimum. For a boat there are six vectors, three 
maximum values and three minimums, These vectors are tested to find 
the positions which are in fact farthest from and closest to the datum 
which served as the starting point for the drift determination. Using 
the maximum and minimum values obtained in the test as entering 
arguments, a table 44 stored in memory produced the final value of drift 


error, d 


a See Figure 10-D. 


The final determination of the total probable error of position, c, 
is then accomplished by combining the three separate errors in accordance 


with the equation: 


c = VV x4 + v2 + a? 


See Figure 10-E. 
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and tne total probable error of the position has been determined, the 
boundary limits of the search area may be defined. These boundaries 
may be expressed in either of two ways. The first is in the form of the 
length of the search radius which extends from the datum point. This 
method is used for vector type searches. The second is in the form of 
geographical coordinates of latitude and longitude defining a rectangular 
search area centered on datum with a length and width each equal to 


twice the search radius. See Figure ll. 
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The specific input data needed to define tne search area is: 

1. The number of the search being conducted. This is needed to 
ascertain what safety factor will be applied in the calculation of the 
radius, R. 

2. The various safety factors associated with each different search 
number. These are the factors listed in Chapter II, page 29. 

3. The total probable error of position, c. This has been computed 
in the preceding phase. 

4, The coordinates of the forecasted datum which were determined 
in the first phase of the overall program. 

The actual computations are fairly simple. The search radius, R, 
is the product of the safety factor and the total probable error of position, 
c. This may be printed out along with the coordinates of the forecasted 
datum to provide the first method of defining the search area. 

In the second method, the latitude increment will be equal to the 
search radius (in nautical miles) divided by 60 nautical miles per degree 
of latitude. This increment of latitude is added to the forecasted datum 
latitude to yield the northern limit and subtracted from the datum latitude 
to provide the southern limit. 

In determining the increment of longitude equal to the search radius, 


the datum latitude was considered to be the middle latitude, Ly, ofa 
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small area Mercator plotting sheet. The equation relating the search 
radius distance in miles to degrees of longitude ig; 49 

Degrees of Longitude = GUsosEy eo 
This increment was then subtracted from the datum longitude to obtain 
the longitude of the eastern limit and was added to the datum longitude 
for the western limit. 

This then completes the basic program for the determination of the 
search area problem. However, the specific program which was developed 
contains three additional items of interest. The first two of these are 
designed to provide general information concerning the possible conditions 
of the survivors. This information assists the RCC Controller in making 
the difficult decision of when to terminate an unsuccessful search and 
rescue case. It will also provide information to the search units 
-indicating tentative survival times. 

The first of these two subroutines interrogates a table stored in memory 
which provides information expressing the effect of wind and air temperatures 
@hH exposed survivors. This table converts the values of wind and air 
temperatures at the scene of the incident to the equivalent temperature 
on exposed flesh if the wind velocity were zero. 46 

4° Bowditch, op. cit,, p. 89. This section describes the various 
mathematical relationships between latitudes, longitudes, and distances 
on small area Mercator plotting sheets. 

49This information was obtained from a "Wind Chill Chart" developed 
by tne U.S. Army. It appeared in the December 18, 1964 Commandant's 
Bulletin, No. 51-64, Supplement No. 4. published by the U.S. Coast Guard. 


This table was reprinted from U.S.N. Medical Newsletter, Vol. 32, No. 12. 
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The inputs needed are two in number. First, the air temperature at 
the scene is needed and may be obtained from weather reports or forecasts 
for the area. The second input is that of the velocity of the local surface 
winds in the area. This has been entered in an earlier phase. 

These two inputs serve = the entering arguments for the table which 


then reads out the desired information directly. See Figure 12. 
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The second oe oweme associated with the condition of the survivors 
deals with immersion hypothermia, thatis, the life expectancy of 
survivors immersed in water at different temperatures. The desired 
information is an estimate of the time survivors will remain conscious 
and an estimate of the expected time they have before death becomes an 


almost certainty. 
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The only input needed for this portion of the program is that of the 
sea water temperature in the area of the distress incident. This information 
is also available in the local weather reports or forecasts for the distress 
area. 

This input information is the entering argument into tables which 
provide the time ranges at different temperatures for consciousness and 
time until death.47 See Figure 13. 

The third and final item pertains to the format of the final output. 

In order to expedite the entire process of transmitting information to the 
searcn units which may be assigned, the format of the computer printed 
output was designed to be that of the standard Outgoing Message format. 
Instead of taking the information from a given computer run, transcribing 
this information to message blanks, and then routing it to the communications 
center, the actual print out serves as the drafted message. The only 
information which must be inserted is the message date-time group (DTG), 
address or name of search unit involved, identification of the distress 
incident, and the releasing officer's approval. This also eliminates the 
possibility of transcription errors in the many numbers involved. 
Appendix B. contains samples of the print out format for the first and 
subsequent searches with computer calculated figures indicated by 
47Climatological and Oceanographic Atlas for Mariners (Washinaton: 
Office of Climatology and Oceanographic Analysis Division, Department of 
Commerce, 1961), Chart No. 157. 
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underscoring and RCC Controller entered information indicated by 
underscored X's (XXX). In the space after the last standard print out 
paragraph the RCC Controller can enter any additional information that 


he deems appropriate. ' 
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CHAPTER IV 


RESULTS 


The specific computer program which resulted from this study is 
found in Appendix B. As each portion of the program was developed, it 
-was tested in such a way that every part of the program was operated. 
Intermediate print-outs of various steps were obtained and verified for 
all possible types of inputs for that particular portion. Efforts were 
directed toward simplifying and shortening the routines involved. 

As the final program was amalgamated, the testing and review 
continued. When it was found to be functioning properly, the various 
.test print-out statements and corresponding comment cards were remov@d. 
Only those comment cards which will assist in reviewing the program ee, 

Sample illustrations of a distress case involving an aircraft (ditching) 
with associated survival craft ranging from life jacket, life raft to life 
boat, with a ship as the search unit are found in the following section. 
These problems show the input data, the manual solution, and the 
computer solution. Manual solutions were done on small area Sain 
sheets. It was noted that when the manual solutions were done by people 
thoroughly familiar with the problem, human error lead to one 17 mile 
error in datum. The computer solution was correct. 


The running time for the samples varied with the number of the search. 
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The longest time, 68 seconds, occurred on the run for the initial search, 
and 54 seconds were required for the subsequent searches. The first 
search has a longer print-out associated with the distress message. 

It also involved computations of initial position error of the distress 
incident and SAR unit navigation error as well as the immersion hypothermia 
and chill factor relationships which were not computed in subsequent runs. 
All rfnning times would be reduced if this program were in an operational 
status as the various tables, program instructions and other associated 
data would already be in storage within the computer. This information 


had to be entered at the beginning of each of these sample runs. 
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I, TRIAL DATA COMPUTATIONS 


RUN ONE 


INPUT DATA 


Classification of distress unit position report 


Classification of SAR unit navigation method 


Initial reported distress position 

Position of distress unit's last fix 

SAR unit position at start of case 

Time of distress incident 

Time specified for arrival of SAR unit 

Surface wind at distress position 

surface wind for past 24 hours at distress 
position 

Average surface current at distress position 

Air temperature at distress position 

Sea water temperature at distress position 

Search number 


Type lifesaving craft employed 


(A 


Aircraft = DR 
Ship = FIX 
42.50N, 52.42W 
41.50N, 61.50W 
50.00N, 35.00W 
301730Z (JAN) 
021130Z (FEB) 


340°T @ 17 kts. 


265°T @ 10 kis. 

195°T @ 15.6 miles/day 
12. 89°C 

(eee 

1 


Life raft 





RESULTS BY MANUAL COMPUTATIONS 


Distress unit navigation error (X) 
SAR unit navigation error (Y) 
Forecast datum of survivors 


Forecast search area limits 


51.3 miles 
5.0 miles 
40.97N, 51.80W 


42.14N, 39.81N 
50. Z5W, Soegonv 


RESULTS BY COMPUTER (PROGRAM) COMPUTATIONS 


Distress unit navigation error (X) 
SAR unit navigation error (Y) 
Forecast datum of survivors 


Forecast search area limits 


RUN TWO 


INPUT DATA 


Last computed distress position 

Type lifesaving craft employed 

Time associated with last distress position 
Time specified for next search to begin 


Surface wind at computed distress position 
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51.1 miles 
5.0 miles 
40.90N, 51.84W 


42.01N, 39. 80N 
90.37W, 53.30W 


40.90N, 51.84W 
Man in water 
0211302 (FEB) 
0413302 (FEB) 


180°T @ 20 kts. 





Surface wind for past 24 hours at computed 
distress position 

Average surface current at computed distress 
position 


Search number 


RESULTS BY MANUAL COMPUTATIONS 


Forecast datum of survivors 


Forecast search area limits 


340°T © 17 kts. 


195°T @ 15.6 miles/day 


Z 


39.88N, 52. 16W 


41.30N, 38. 46N 
00.28W, 54.04W 


RESULTS BY COMPUTER (PROGRAM) COMPUTATIONS 


Forecast datum of survivors 


Forecast search area limits 


RUN THREE 


INPUT DATA 


Last computed distress position 
Type lifesaving craft employed 
Time associated with last distress position 


Time specified for next search to begin 
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39.88N, 52. 14W 


41.30N, 38. 47N 
50.30W, 53.98W 


39.88N, 52. 14W 
Life boat 
041330Z (FEB) 


0621302 (FEB) 





Surface wind at computed distress position 095°T @ 04 kts. 
Surface wind for past 24 hours at computed 

distress position 180°T @ 20 kts. 
Average surface current at computed distress 

position 180°T @ 15.6 miles/day 


Search number 3 


RESULTS BY MANUAL COMPUTATIONS 

Forecast datum of survivors 39.89N, 52.09W 

Forecast search area limits 41.62N, 38.16N 
49.83W, 54.35W 

RESULTS BY COMPUTER (PROGRAM) COMPUTATIONS 

Forecast datum of survivors 39.87N, 52.07W 


Forecast search area limits 41.60N, 38.14N 
49.82W, 54.32W 


SUMMARY OF RUNNING TIMES 


Run number one l min. 8 sec. 
Run number two QO min. 54 sec. 
Run number three 0 min. 54 sec, 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 


I, SUMMARY 


Growing maritime operations have found the need for more effective 
search and rescue facilities expanding. One area of search and rescue 
operations, that of the determination of search areas and datum points, 
has been studied to ascertain the feasihility of the use of computers in 
solving this problem. Using the National Search and Rescue Manual as 
a guide to the present doctrine employed by this country, a computer 
program was developed, tested and found to be an operationally feasible 
program. 

During the course of the development of this program certain areas 
of the current doctrine were reviewed in the light of recent literature and 


theories which may offer improved techniques. The conclusions reached 


in this review are discussed in the following section. 
Il. CONCLUSIONS 


The basic doctrine outlined in the National Search and Rescue Manual 


is amenable to computer programming. Although the program which was 


developed does not cover every possible aspect of this doctrine, it does 
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handle one of the basictypes of search problems frequently encountered. 
This program can greatly reduce the time required to solve the problem 
associated with the location of datum and its surrounding search area. 
The accuracy of the computer solution is dependent entirely upon the 
accuracy of the input information. It greatly reduces the possibilities of . 
human error in carrying out the rather extensive calculations involved in 
the problem, As new methods are developed for determining various 
portions of the input data with greater accuracy, the value of the computer 
solution will greatly increase. 
conclusions were reached. The first of these was that the doctrine 
outlined, although not the ultimate in sophistication, certainly represented 
an excellent model for use in existing RCC facilities which were not 
equipped with automatic data processing equipment or with large staffs. 
Admittedly there were areas of weakness in the basic doctrine. However, 
to refine the doctrine to fully consider the different shortcomings would 
have unnecessarily complicated the problems at the sacrifice of valuable 
time. And there is no definite measure of how this additional refinment 
would alter the risks associated with the resultant decisions. 

However, the review of the literature in this study did suggest 
several possible areas for improvement if a computer installation is 


available, and one suggestion is applicable with existing manual operatéons. 
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Considering the computer applications first, it is apparent that the 
problem of free fall and parachute drift is subject to operational analysis 
by computers. The present running time of slightly less than ten minutes 
associated with Cohan's program can probably be reduced as his program 
is refined. 

Also, a computer could use the program developed by Burns at the 
Navy Oceanographic Office to determine a drift elipse area based upon 


the statistical current roses found in the Atlas of Surface Currents. 





The development by the Fleet Numerical Weather Facility of the 
system for forecasting surface currents appears to have the greatest 
potential for improving present doctrine in that it may be used either with 
or without a computer. The program developed by this study can be 
modified to receive information from the Fleet Numerical Weather Facility 
on.the predicted surface currents at the four points of the grid containing 
datum. These vectors, which reflect a combination of the permanent 
circulation pattern and local wind currents, can then be resolved into a 
vector for the datum point. Since this information is based on recently 
reported weather conditions and sea surface temperatures, it should tend 
to be more accurate than the information taken from the Atlas of Surface 
Currents which represents a statistical summary of many years based, in 
some cases, ona single observation. 

This analysis of Fleet Numerical Weather Facility information cquld 
also be graphically determined or manually calculated without the assistance 
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of a computer. 

Laevastu's work indicates that local wind current divergence can be 
treated in a more refined manner by considering small increments of 
latitude and wind velocities. This refinement could be handled very 
easily by modifying that portion of the specific program relating to this 


computation. 


Ill, IMPLICATIONS 


The program developed in this study can be rewritten in the computer 
language used in the AMVER Center, Office of the Commander Eastern 
Area. Thus, it can serve to assist the handling of these distress cases 
which occur in the North Atlantic. It is suggested that serious consideration 
be given to the adaptation of the Fleet Numerical Weather Facility data 
services and the associated modification of this basic program. Should 
this latter suggestion not be practical, consideration should be directed 
toward utilizing Burns' program in conmnenion with the Atlas of Surface 
Currents, As the Pacific AMVER system becomes fully operational, this 
this program can be further modified to handle SAR cases falling under the 
jurisdiction of Commander Western Area. 


In view of the number of entries found in the Atlas of Surface Currents 





(approximately 20, 000 for the Atlantic and 40, 000 for the Pacific) it is not 


recommended that this information be placed in storage and the basic 
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program modified to handle that portion of the Average Sea Current problem 
automatically. Individual cases should continue to use this information as 


manually extracted from the appropriate charts. 
IV. RECOMMENDATION OF AREAS FOR FURTHER STUDY 


This study has served as a first step in the development of a full-scale 
computerized SAR program. As such, it has opened up many new areas where 
further study and research should prove fruitful. Some of these areas are 


discussed briefly below. 


incident and estimated position are not known. 





When a ship or aircraft is missing or overdue, the RCC Controller 
is faced with the problem of estimating the probable boundaries of an 
area in which the distress may have sae He must establish the 
limits of the time period in which the distress could have occurred. With 
this informationhe:then.applies minimum and maximum possible values of 


drift for the time period in question starting with the earliest possible 


distress time and ending with the latest possible distress time. 


Development of a fully operational free fall and parachute drift program. 





Refinement and further development of the work done by Armacost, 
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Saunders and Cohan can materially contribute to the effectiveness of a 


computer program for SAR operations. 











Development of a program to handle small craft (inshore) distress cases. 


With the increasing interest in recreational boating displayed by 
the general public, it is mandatory that improved techniques be developed 
to handle the associated increase in SAR cases occurring with these craft. 
It has been estimated that approximately 85% of all SAR cases handled by 
the U. S. Coast Guard occur within twenty to twenty-five miles of the 
coast. Someplace in the neighborhood of 95% of all cases are found to 
take place within 100 miles of the coastline. 48 

One of the major difficulties in treating these cases is associated 
with the effects of tidal currents, rotary currents, and run-off currents, | 


all of which are effected by the configuration of the coastline. 


Development of a program relating search area, search units, and search 





patterns. 


Once a search area has been specified, the ROC Controller is faced 
with the problem of recommending search patterns for those available 
48 R. House, Analysis of U. S. Coast Guard Aviation Operations 


(Report Number SER=50371, Stratford: Sikorsky Aircraft Division of United 
Aircraft Corporation, 1964), p. 36. 
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search units which are participating in the case. Since the type of search 
pattern recommended is a function of the number and type of search units 
participating, the weather and other operating conditions at the scene, 
and the various probabilities of detection desired, the development of a 
program to assist in making this recommendation would be most beneficial. 
In essence this would be a program treating the doctrine outlined in 


Chapter Seven of the National Search and Rescue Manual. 


eat = 2 feet Se SE Ee SS eee oe 
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facilities. 


As the new system of reports for distress cases becomes fully 
operational, actual historical data will become available. Analysis of 
this data and existing SAR facilities should shed new light on possible 
relocation and expansion of SAR facilities. Studies similar to that performed 


by DuPeza, Eustis, Guthrie and Zook49 should be of great value. 


Continued drift analysis. 





At the present time the U. S. Coast Guard has a drift datum marker 


boy which is being evaluated by operational field testing. The data 


4 
‘ 


-obtained from these tests coupled with the current studies being conducted 


43 tules B, DuPeza and others, "A Method for Evaluating SAR Requirements 
for the Twelfth Coast Guard District" (Unpublished Master's research paper, 
The United States Naval Postgraduate School, Monterey, 1964). 
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by the U. S. Naval Oceanographic Office, U. S. Naval Fleet Numerical 
Weather Facility, and other oceanographers should provide a great deal 
of information which can lead to a better understanding of the problems 
associated with ocean current analysis. 

In all cases, the continued application of man's growing knowledge 
and utilization of the tools which he has developed can be directed 
toward the benefit of mankind. The success or failure of such efforts 
is limited only by the imagination, enthusiasm and patience of those 


who truly seek the answers. 
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APPENDIX A. 


FLOWCHART SYMBOLS AND DEFINITIONS 
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FLOWCHART SYMBOLS AND DEFINITIONS 


Terminal indicator 


- Input/output indicator 


pe: - Decision indicator 


- Processing, annotating indicator 


- Off page connector 


C) 2 & On page connector 


== fy - Flow indicator 
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APPENDIX B. 


PROGRAM SAR 


Program SAR is the specific program developed in this study. It 
determines the datum point and associated search areas for distress 
cases involving ships and aircraft which ditch. This Appendix includes 
the definitions of the terms and inputs found in the program, the 
associated computer program itself, and sample print-out as discussed 


at the end of Chapter ITI, 


on 


VARIABLE 


NAME 


ATEM 


CRAFT 


DCLS 


DDATE 
DFTERR 
DHR 
DLAT 
DLONG 
DMIN 


DMO 


DEFINITIONS OF ALL INPUT VARIABLES 


qaqa ae le ee ee 


DEFINITION 


Air temperature at distressed position in degrees centigrade 
(CGH 


Type lifesaving craft that ROC Controller thinks survivors 
will employ: 
Life jacket = code l. 


Boat - code 2. 
Raft - code 3. 

Classification of Distressed unit's position report: 
Surface DR = code 01 
Aircraft DR - code 02 
Surface FIX = code 05 
Aircraft FIX = code 10 
Surface and Aircraft DF 

ClassA - code 05 
Glass B = code 20 
Class © - code 50 
Surface and Aircraft DF-FCC 
Class A -~ code 20 
Class B - code 40 
Class C - code 60 
Class D = code 80 


Date of distress in ZULU (GMT) time. 

Drift error from Min-Max Table in CG-308. 

Hour of distress in ZULU (GMT) time. 

Latitude of distress in degrees and hundredths. 
Longitude of distress in degrees and hundredths. 
Minute of distress in ZULU (GMT) time. 


Number of days in the month associated with DDATE. 
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DROG 


FLAT 


PWDIR 
PWIND 


SCLS 


SCURD 
SCURV 
S DATE 


SEARCH 


Drift of a raft or boat without a drogue as a factor of wind 
velocity. Values obtained from a graph in CG-308. 


Time until an immersed man will become exhausted or 
unconscious as a factor of sea water temperature. Obtained 


from table in Atlas of Surface Currents, 


Latitude of last fix of distressed unit in degrees and 
hundredths. 


Number of the search that is being computed. 


Wind direction at the distress position for the 24 hour period 
prior to the distress. 


Wind velocity at the distress position for the 24 hour period 
prior to the distress. 


Classification of the SAR unit's position determination: 


Surface DR - code 03 
Aircraft DR - code 
Surface FIX - code 05 
Aircraft FIX = code 10 
Surface and Aircraft DF 
ClassA = code 05 
Class B = code 20 
Class © = code 50 
Surface and Aircraft DF-FCC 
Class A - code 20 
Class B - code 40 
Class © - code 60 
Class D ~- code 80 


Direction of surface current at distress position taken from 
the Atlas of Surface Qurrents. 


Velocity of surface current at distress position taken from 
the Atlas of Surface Gurrents. 


Date specified by RCC Controller, usually the ETA of the SAR 
unit on scene in ZULU (GMT) time. 


Safety factor associated with 'M'. 
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SHR 


SLAT 


SLOW 


SMIN 


STEMP 


STM 


W DIR 


WIND 


WINTEM 


X 


Y 


Hour specified by RCC Controller, usually the ETA of the SAR 
unit on scene in ZULU (GMT) time. 


Latitude of the SAR unit's last fix in degrees and hundredths. 
Longitude of the SAR unit's last fix in degrees and hundredths. 


Minute specified by the RCC Controller, usually the ETA of 
the SAR unit on scene. 


Temperature of the seawater at the distress position, OC. 


Survival time for a man immersed in water at a certain 
temperature as a factor of the temperature. 


Wind direction at distress position in degrees True. 
Velocity of wind at distress position. 

Wind chill factor for the effect of wind on exposed flesh. 
Position error for the distress report. 


Position error for the SAR unit. 
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VARIABLE 
NAME 
“DCLS 
SCLS 
DLAT 


DLONG 


PLAT 


FLON 


SO LAT 


SLON 


CRAFT 


DDATE 


DHR 


INPUT DATA FOR FIRST SEARCH 


DATA CARD 
COLUMNS 


Lae 


ll=16 


WAT 


Zoaeo 


29-34 


35-40 


41, 42 


43, 44 


45, 46 


SOURCE AND FORMAT 
OF INPUT DATA 


Classification of distressed unit's position 
report. Two digit code number. 


Classification of search unit's position report. 
Two digit code number. 


Latitude of distress incident. Six digit number 
of degrees and hundredths (i.e. 45939' = 045.65). 


Longitude of distress incident. Six digit number 
of degrees and hundredths (i..e. 75921'= 075.35). 


Latitude of last fix of distressed unit six digit 
number of degrees and hundredths. (i.e. 
42°95 4'= 042.90). 


Longitude of last fix of distressed unit. Six 
digit number of degrees and hundredths (i.e. 
63°13'= 063. 22). 


Latitude of SAR unit at start of incident. Six 
digit number in degrees and hundredths (i.e. 
40°38'= 040. 63). 


Longitude of SAR unit at start of incident. Six 
digit number in degrees and hundredths (i.e. 
68°07'= 068.12). 


Type of life saving device that the ROC coordinator 
believes survivors will be using. One digit code 
number followed by a decimal. 

Date of distress incident. Two digit number from 
the distress DTG. 

Hour of the distress incident. Two digit number 
from the distress DTG. 
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VARIABLE 


NAME 


DMIN 


S DATE 


SHR 


S MIN 


WIND 


PWIND 


WDIR 


PWDIR 


SCURV 


SCURD 


DATA CARD 
COLUMNS 


47, 48 


49,50 


91,52 


93,94 


99, 96 


07,98 


99-61 


62-64 


65-68 


69-71 


SOURCE AND FORMAT 
OF INPUT DATA 


Minute of the distress incident. Two digit 
number from the distress DTG. 


Date specified by the RCC Controller, usually 
the date (GMT) that the SAR unit will arrive on 
scene. Two digit number. 


Hour specified by the RCC Controller, usually 
the hour (GMT) that the SAR unit will arrive on 
scene. Two digit number. 


Minute specified by the RCC Controller, usually 
the minute that the SAR unit will arrive on scene. 
Two digit number. 


Wind velocity at distress position from the latest 
weather information. Two digit number. 


Wind velocity for the previous 24 hours at the 
distress position. Two digit number. 


Direction from which the wind is blowing at the 
distress position at time of distress. Three 
digit number. 


Direction from which the wind had been blowing 
at the distress position for the 24 hour period 
prior to the distress time. Three digit number. 


Surface current velocity at the distress position 
taken from the Atlas of Surface Currents. Four 
digit number (nautical miles per day) in miles 
and tenths (i.e. 09.9). 


Direction of the surface current flow at the 
distress position taken from the Atlas of Surface 


Currents. Three digit number. 





VARIABLE 
NAME 


ATEM 


STEMP 


DMO 


DATA CARD 
COLUMNS 


12275 


76-79 


80 


SOURCE AND FORMAT 
OF INPUT DATA 


Air temperature at the distressed position at the 
distress time taken from weather reports in degrees 
centigrade. Four digit number in degrees and 
temtusmine. 1746). 


Sea water temperature at distress position in 
degrees centigrade taken from latest weather 
data. Four digit number in degrees and tenths 
(ijn 1 Sana). 


Search number for which datum is to be computed. 
One digit number. 


Number of days in the month of DDATE, Two 
digit number. 


Left blank for the first search. 


Left blank for first search. 
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VARIABLE 


NAME 
DCLS 
SCLS 


DLAT 


DLONG 


FLAT 
FLON 
SLAT 
SLON 
CRAFT 
DDATE 
DHR 
DMIN 


S DATE 


SHR 


SMIN 


DATA CARD 
COLUMNS 


vee 
3,4 


Sa llig 


11-16 


IAPs 
Z3—20 
29-34 
35-40 
41,42 
43, 44 
45-46 
47, 48 


49,50 


oro 


03, 54 


INPUT DATA FOR SECOND THRU FIFTH SEARCHES 


SOURCE AND FORMAT 
OF INPUT DATA 


Blank. 
Blank. 


DTMLAT from preceding search results. Six 
digit number in same format as DTMLAT. 


DTMLONG from preceding search results. Six 
digit number in same format as DITMLONG. 


Blank. 

Blank. 

Blank. 

Blank. 

Same as CRAFT in preceding search. 

SDATE from preceding search. 

SHR from preceding search. 

SMIN from preceding search. 

New specified date by RCC Controller, usually 

the estimated date (GMT) that the SAR unit will 
complete the preceding search. 

New specified hour by the RCC Controller, usually 
the estimated hour (GMT) that the SAR unit will 
complete the preceding search. 

New specified minute by the ROC Controller, usually 


the estimated minute (GMT)that the SAR unit will 
complete the preceding search. 


98 





VARIABLE DATA CARD SOURCE AND FORMAT 


NAME COLUMNS OF INPUT DATA 
WIND 55, 96 Wind velocity at the preceding computed datum 


position. Taken from weather reports for the 
preceding specified time (SDATE, SHR, SMIN). 
Two digits in same format as search one. 


PWIND 97, 58 Wind velocity for the previous 24 hours period 
at datum computed in the preceding run. Two 
digit number obtained from weather reports. 


WDIR 59-61 Wind direction at datum computed in the preceding 
run. Three digit number obtained from the latest 
weather report. 


PWDIR 62-64 Direction of previous 24 hour wind at datum 
computed in preceding run. Three digit number 
obtained from weather reports. 


SCURV 65-68 Surface current velocity at datum computed in 
preceding run. Four digit number of miles and 
tenths (i.e. 13.4) obtained from the Atlas of 
Surface Currents. 





SCURD 69-71 Direction of surface current at datum computed in 
preceding run. Three digit number obtained from 
Atlas of Surface Currents. 


. ATEM Toews Blank. 
STEMP 76-79 Blank. 
M 80 Number of the search, 2 thru 5. One digit. 
DMO LZ Number of days in the month of DDATE. Two 


digit number. 


X 3-7 Value of 'X' from the initial search. Five digit 
number with one decimal place. 


Ye 8-12 Value of'Y' from the initial search. Five digit 
number with one decimal place. 
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FIRST SEARCH PRINT OUT 





OO 060,0 Ys 


FM COMEASTAREA 


TO XXXXXXXXXXXXXXXKXXXKKK 


Bil 


UNCLASS 


A. DISTRESSED XXXXXXXXXXXXKXX 


ly 


2 


FORECAST DATUM AT 021130Z. 40.90N, 51.84W. 
FORECAST SEARCH AREA BOUNDED BY 42.01N, 39.80N, 
50.37W, 53.30W. 

SURVIVOR CONDITION FROM 3017302 OF MAN IN WATER. 
TIME UNTIL EXHAUSTION 1.00 HOURS. TIME UNTIL 
DEATH 3.00 HOURS. 

EFFECT OF SURFACE WIND SAMEAS -27.0 DEGREES 


FAHRENHEIT WITH ZERO WIND. 


RELEASED BY XXXXXXXXXXXXXXXXKXKK 
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SUBSEQUENT SEARCH PRINT OUT 
Xeoeese. 47 
FM COMEASTAREA 
TO SOOOOOREEOSUCOCOOOO 
BT 
UNCLASS 
A, DISTRESSED XXXXXXXXXXXXXXXXX 


1, FORECAST DATUM AT 0621302. 39.87N, 52.07W. 





2. FORECAST SEARCH AREA BOUNDED BY 41.60N, 38, 14N, 





RELEASED BY XXXXXXXXXXXXXKKKXKXXXX 


101 





STl 


vIT 


etl 


cll 


Ee 


OTT 


60T 


80T 


EOE 


90T 


SOT 


7OT 


COT 


cOT 


TOT 


& 


(L°T=1*(1)WiS*(1)HX9) £006 aVay 
(L)WLS*(L)HX NOISNAWIG 

(T°94TT) LVWYOS 

(HTS THIS (TTS Tere (rst) WSLNIM)) *S6T avay 
(TTS T)WALNIM NOISNAWIG 

(T*€4S) LVWHOS 

($*T=4*(¥)HDYVaS) “TOs avay 

(S)HOYVIS NOISNAWIG 

(0*Z46T) Ly¥WOos 

(GTST=14 (GTS Lar e( rs 1) uyaisa))*2tet avay 
(ST*ST)YYILIG NOISN]WIG 

(T*h4GTt /T%748T) LVWYOS 
(€€*T=N¢(N)9ONG) *2T% avaY 

(€€)90NG NOISNAWIG 


uVS wWVY90dd 


Gl 


ot 


To08 
008 
EOS 
eal 
TLV 
O) taal 
liv 
tty 


OV 


aey, 





A! 
a ae | 
Len 
Oel 
6cl 
8c 
Let 
921 
Gcl 
721 
Ect 
CEN 
UGE 
O2T 
6TT 
8tT 
LTT 


ITT 


NO1IS=NO1Z 
i IVIS=1V 12 
SHs-5 1D C0 | 
VOTSEOTSEOT(*%4¥-S1DS)4I SOT 
S1d-=X eeot 
L6T OL O9 
NO1V4S=NO1Z 
1V14=1V12 
Sum d=S 1) TOT 
ZOTS TOTS LOT(*2-S1Dd0) JI Oot 
*(4) YOUN NOTLVOLAWN LINN UVS GNV (X) YONYNS NOILVOIAVN LINN 
SSIYLSIG 4O SANIVA SHL JLNdWOD TT HONOYHL OOT YSEWNN SLN3SWSZLVILS 
€07*O0TSOOT(I-W)4I 
(1°G42*0°Z4/T1 Led *0* Ed L 4460*E4260°24662°94960°Z4Z) LYNNOS LOOE 
ASX*SOWGSWS dWILSZ 
SWALVSAUNDS AYNIS *YLAMd *YIGMSGNIMd*GNIM*NINS*YHS* 3LVGS *NIWGSUHGT 
‘3lvads LAVYD*NO1S* LVIS*NO14*1LV14*9ONO1G* LV10*S19DS*S1Dd* TOOE AVAYOOOOE 


esse eo aia bb) | VWNOS 006 


i f0¥6 





tst 
OST 
67T 
871 
Lyt 
a7 
GT 
71 
ert 
cyt 
Tv 
O7T 
6¢eTf 
SET 
Lol 
rae || 
Seu 


Uae | 


(3INVL)4NVLV=aD 
W4IIG/O1410=32NVL 
(WLVIZ-WLV10) 4Sav=WsI 
(°O8T/9THT°ExLVIZ) INIS#699Z*°E2-TZNVISeL *ST6L=WLVIZ 
(°OST/9THIT* Ex L VIG) INIS%689Z2°EZ2-TONVL4¥L *ST6L=WLYIG 
(TZNvL)4907%6246%°O=1ZNVL4 
(TGNVL)3907%627€%°0=10NV 14 
( (1ZONV) 4S0D/ (1ZONV) SNIS) 4SQV=7ZNVL 
( (JGONV) 4S0D/ (1G9NV) SNIS) SSAV=70NVL 
-OST/OT HIE %( °Z/LVIZ+° SH) =1Z9ONV 
"ORT/9THT° Ext *2/LV 104° SY) =109NV 
( LV1Z-LV10) 4SaVx°09=V1410 
(NO1Z-9NO10) SSGV¥*09=015I0 L6T 
“SNOTLVLNdWOD SLUVd TVNOIGIUAW AG 
SNOILISOd OML WOU4 ADNVLSIG JLNdWOD LOT OL L6T YASWAN SLNAWALVLS 
| E04 OL 09 
SloS=A 7OL 


EGE SG .0O9 


B 


104 





69T 


SOT 


L9T 


pee) | 


GOT 


79T 


cot 


Cot 


ton 


O9T 


6Gl 


Sot 


LGt 


9GT 


GGT 


9ST 


eo 


cSt 


S07 eOn OS 


SYUHS+SUYHG—- ° 2% (ALVGS+3LVGG-OWd) =aWIL 


T0L*G07'G047(341LVOS—-31LVdd) JI 


°“O9/NIWS+YHS=SUYHS 


°“O9/NIWO+UHOG=SYHO 


“UV SHV “a oaeCda eel st ONO Cae oso Ss ldrant 


LWHL SYNOH JO YSEGWAN SHL JLNdWOD SO” HONOYHL €O% YAGWAN SLN3WSLVIS 


(J2)35097 yids °+°Ot=A 
SOOT US 

(39) 43S0D/V1410xG0°+°SG=A 
GOT OL O9 

(39) 43SO0O/V IS TOxT®+°OT=X 
sOl-oOr OS 

(39) 5S09/V 13 10%S0°+°S=X 
Onaga Wl eee oo) 4) 


7 oe LL ee el) ) J] 


So alonttse 


SSO SOs 0-515) 3! 


ZOL 
ToL 
OoL 
70% 


COV 


7TT 


et 


et I 


Pet 


OL 


60T 


gO 


LOT 


109 





e 
La a I al 
« 


LST 


98T 


GBt 


78T 


est 


CRT 


T8l 


OST 


6LT 


Sil 


LLT 


ILT 


GdT 


7LT 


Cra 


CA 


tet 


OLT 


GVUM—"Z/9THT° Sx? G=lGvy 
OST/9THT* ExdIGM=aVvUn 

°O4-YIGM=YIGM 

€*T=N 124 Od 

°O9Z+NIGM=YIGM 

(EINIWAS(EXWWAS (EINTNX€ (€ 1 XVWX NOISNI3WId 
°42/(GNIMN)90Nd=33 10M 
€€=CNIMN 

00046 hI % yLy(€E-GNIMN) JI 

GNIM=GNIMN 

°72/(L8°T/(G°2+GNIM) )=34910 

Ely OL 09 

*¥2/ONIMx8°=3310 

60468076 80%(°S-GNIM) SI 

*AVM3391 4O 

SLN3NOdWOD A GNV X JHL JLNdWOD LZ% HONOYHL LO” YISWAN SLNAWIFLVILS 
LOv S108 82H °2-L4VYD) 4! 


Sato ocmGmeic + Ca byOdd]4iVdS ))=3aWI 1 


OV 
6Tv 
8% 
Lt? 
ote 
GT” 
VT 


000% 


ely 


6047 


80% 


LOVv 


90% 


GOV 


106 





ee eee eee 


TA ts ete ee 5 


GO? 
70e 
cO02¢ 
COZ? 
T0¢ 
OO 
66T 
8é6T 
Lod 
961 
G6T 
76T 
cét 
col 
Tél 
O6T 
6st 


8st 


AWI Le (°O8T/9TYT*ExXGYNIM) ~* Z/9THT °E%°S) INISKUNDMEYNDA 
“BWI Le( (*O8T/9THT°€xGYNIM) -°Z/9THT PE %°S) ISOD¥YNDIM=UNIX 
O€42 OL O09 

°OLZ+YIGMd=GyNoM 

O€4%2 OL O9 

°09€-GYNDM=aUNNoM 

LEZ * Een SEE ( °09E-GYNDM) JI 

-O8T+YIGMd=ay"noM 

c€7 SOE 0€4( °OT-LVI0) JI! 

°yZ2/(8T*T/GNIMd) =yNoM 

: 14140 LNAYYND GNIM WvD07 sO 
SLNANOdWOD A GNV X JHL 31NdWOD bE HONOYHL 8247 YIGWAN SLNAWALVIS 
“Ov T-YIGM=YIGM 

ANNI LNOD 

3WI Lx (1 dVu) SSOD%3aTOM= (N) XVWX 


JWT Lx (1Td0V8) sSOD%443 10=(N)NIWX 


IWITLx(TdQVY) SNISx34dT0M=(N)XVWA 


IWITLx(TOQVY)ANISx¥3310=(N)NIWA 


HEY 
EE 
ZEbT 
Zen 
ZEHZ 
Lenz 
O€72 
O€h 
62h 


8¢% 


LGAs 
z£acV 
UE) 
Ccy 
Gou 


ay 


107 





a 


ee tn me 


C2 


COC 


PAA 


O22 


602 


‘Sid 


r< 


9T2 


Ste 


aT? 


EA? 


Ce 


tid 


OT 2 


602 


802 


LOZ 


902 


(N)XVWLAG=(N)LATYC 

(2% (N)NIWYOA+Z%%(N) NIWYOX) SLUOS= (NINIWLAO 
(Zxx%(N)XVWYHGA+Z%%(N) XVWYOX) SLYOS=(N) XVWLAG 
(N)NIWA+YNDA+YNDSA=(N)NIWUGA 
(N)XVWAFtYNDAFYNDSA=(N) XVWUGA 
(N)NIWX+YNDX+YNDSX= (N)NIWYGX 

c | (N) XVWX+YNDX+YUNDSX=(N) XVWUAGX 
E*T=N G9” OG 

° 141Yd 40 SLN3NOdWOD WAWINIW 


GQNV WOWIXVW TLOL ALNdWOD OL% HONOYHL 8S%7 YSEWNN SLNAWSLVIS 


(EC YNIWYGQA © (E)XVWYHGAT 


79% 
C77 
Zo 
how 
09% 
6S% 


8S47 


“(CEYNTWYHCXS (CE) XVWYHAXs (CINIWLAGS (€) XVWLAGS (9) LAST NOISNAWIGOLSY 


JWI Le ( *O8T/9T7 T° €xGdNDS—-7S8°L) SSOD%* 72/AYNDS=UNDSX 
JWI Lx (°O8T/9OTHT PE xOYNDS-7S8°L) INIS%°72/AYNDS=YNISA 
Todos 

4G%61G6HSEG7(AYUNDS) AI 

°L4ITYd LNAYYND JDVAYNS JOVYSAV AO 


SLNSNOdWOD A QNV X 3HL 3SLAdGWOD GS HONOYHL 267 YSGWNN SLNSWSLVILS 


GSGv 


7G 


tSyv 


CG” 








—— EEE 


T7e 


Ove 


6E¢ 


BEC 


Led 


Jee 


GE? 


Vee 


A 


Coe 


Lee 


O€2 


62¢ 


Sac 


Lec 


KEG 


Gc¢ 


Ge 


6976976024 ((N) LITYO-NWLSG) JI 
| 9*°Z=N OLY OG 

(T) LIIYG=NWL 4a 

JANI LNOD 

(N) LITYG=xWLsa 
894768946194 ((N) LSIYG—XWL4d) JI 
9*Z=N 894% Od 

(T) L41yG=XW1Jd 

vi% OL 09 

(GS) LJTYG=NW1 4a 

(Z) LITYG=xwWLlsa 

bLY OL 09 

(Z) LITYG=NWL4a 

(G) LdTYd=xXwWLlsa 
997€*997E69907((G)LSTYG-(Z)1L4514NG) SI 
99416997699 ( *2-1L4VYD)AI 

3NNILNOD 


(NINIWLIG=(E+N) LAITY 


89” 


no? 


29% 


S97 


Bore 


eiehye! 


G9v 


109 





6S¢ 


8S¢ 


LG¢ 


9S¢ 


GG¢ 


7G¢ 


ESC 


Coe 


TS¢ 


OS2 


672 


872 


Ly2 


94¢ 


GC 


772 


ed 


Cre 


“OST/9OTHT°€x(0*2/1V10+0°S7) =10ONV 

| *wnivd 

4O JGNLIONOT GALDSFOYNd ALNdWOD STS HONOYHL 7OS YASWNN SLNAWSLVIS 
(*O9/WNLVGA) FLV I0=LVIWLG 

(2x WNLYGA+Z¥xWNLVOX) JLYOS=WLGLSIG 
(ZYNIWYGA+* 7/0 (0 2)NIWYGA—(2) XVWYGA) =WNLVGA 

(2 NIWYGX+*2Z/( (ZINTWYAX= (2) X¥VWYAX) =WAL VAX 

; *WNLVG 

4O 3GNLILV7 G3LDaCONd ALNGWOD €0S HONOYHL 006 YAGWAN SLNAWALYILS 
(c*1)Yyyslsq=sa 

| pap 
SL¥*8LH2s8lH72 (fF)! 
*OoT/NWL3d =f 
"OT/(°O0T+XWL4d) =! 
*(Jqg) yous Ldlyd 
XYWINIW JO 3NIVA JHL JLNGWOD BL» HONOYHL 4¥L YAGWAN SLNAWILVLS 
JANI LNOD 


(N) LITYG=NW1LJa 


70% 


COs 
OS 
Tos 


00% 


110 


Sly 
8172 
LLy 
914 


ey 


OL” 


697 





Or ae mgt th a Ge IE AR le 


nero 


Lie 


Jise 


Gle 


ise 


cic 


CLC 


Lie 


OL? 


ove 


89 


L9¢ 


99°C 


G9? 


09¢ 


Dee Ke 


C7¢ 


Te 


O9¢ 


*SLIWIT V3YV HDYVIS 3LNdWOD L109 HONOYHL €09 YSEwAN SLN3WSLVLS 

| D*¥(W)HDYVSS=4 

°(Y) SNIGVY HDYVIS SSLNdWOD 209 YSGWNAN LN3ZWILVLS 

(2¥x 504+ ZxxA+Z%%X) SLYOS=D 

*(>) NOILISOd 3O YOUNS 31IGVGOUd AVLOL SALNdWOD 009 YSGWNAN LNAWSLVILS 
°09/010WLG-ONO1d= SNOAWWLG 

/ 009 OL 09 

9NO1G=9NO1W LA 

009 OL 09 

°09/010WLG+9ONO1G= SNOIWLG 

GIS *VTSSETSIWNLVGX) JI 

(WALVGA/S (A 2Hx (WLOLSIG/WALVGA)~*°T)4LYOS) ) ISGVHWLGLSIOxWdOSIG=0710W Ld 
(WIWLG-WLV10) 4SavV=ewa4ia 

( *O8T/9TY T° Ex LVIWILG)ANIS%689Z°E2—-WWLONVL¥L°ST6L=WIWLG 

(°OST/9OTHI° Ex LVIG) ANIS%6892°E2—-TONVL¥L°ST6L=WLV 10 

( ( (JWLOONY) SSOD/ (WLGONV) ANIS) SSGV) 4901%627E7° O=IWLONVI 

(( (IGONV) SSOD/(IGONV) ANIS) SSGV) 4901%627E%* O=1GNVL 


7081/91 7 Wee (057 TV IWLGFO.S7)= IWLOONV 


tO he) 


009 


SS 


71S 


ETS 
Cus 
TTS 
OTS 
60S 
80S 
LOS 
90S 


GOS 


HALAL 





nee 


76¢ *VIWYSHLOdAH NOISYsSWWI SALNdWOD 906 HONOYHL OS YSEWNN SLIN3AWSLVIS 


} 


€62 \ ANNI LNOD 
Z6Z 09°04¢8H( (SIT WILNIM-W3LYO4)) SI 
16Z bT*Z=1I 84% Od 
062 °ZE+°G/ *6xWALV=WILYOSD 
682 tt=r 
gaz ‘ ANNI LNOD 
182 0€*0%7607( (7 *T)WSLNIM-NIMYO4D) JI 
982° Il*t=r O€ 0a 
G8Z Ze $ZE $8Z(00°SY-NIMYOSD) JI 
482 GNIMKZET® T=NIMYOSD 
€8Z °YOLDV4 TIIHD GNIM JHL JLNdWOD 8% HONONHL 92 YSEWAN SLN3WILVLS 
282 g09*92*9Z(T-W) JI 
18Z \NO14SG+9NOWLG=9NO1M 
108Z \NO140-9NOWLG=9NO713 
08z °09/( (°O8T/9TYT*ExLVIWLG) 4SOD/4Y) =ONO14G 
612 °09/Y-LVIWLd=LV1S 
Bure °O9/Y+LVIWLG=LV INV 


84 
447 
a 
04 
ZE 
O€ 
62 
gz 
jeg 


2¢ 


Qe: 
707 
SO9 
709 


coo 


Tle 





Cle 
Die 
Oe 
60¢€ 
80€ 
LO0e 
90€ 
GO¢ 
7O€ 
COL 
cOt 
TO¢ 
O0¢ 
66¢ 
862 
EOC 
962 


G6¢ 


JONNOG VANV HDYVSS LSVIDYOS *ZHEE&X66//°MHZ6Z* 146 NHZEZ*L 46 *ZHZ69LZ 
° IV WNLVG LSVOYOS *THOZ*X6*//*G3SS3YLSIC WHET *®X9*//*SSVIDNNHZL & X91 
JIPAGACSXOSSSPOLHCEXGS 7 / OVA SOD Waly le 7 2 lee ee ev td SOG 
ee eNO SIN LOS LY AW LG (ie) Sete eS 1 
(0 )RMS* 9) lOr* ONO IM* ONO Ta fv is i VaINy ONO IWIOS 17 IW1LG Sida ecG INIdd 206 
B06 wo e206 tL L-Wos | 
NIWO+(*O0T#YHG)+(°0000Tx31vVdG)=910r 609 
NIWS + ( ° OO lsdnS) 4 ( OCC Le atvdS J=910A 809 
A-8=A 906 
AJNNILNOD 706 
°OT-~dW3l=dwWw3l S06 
SOG OUG eae ocd Ji aguiails ) 31 
o> l=35 706 Od 
°O8=dW4dl £06 
£06 OL O9 
T=xA 206 
BOG Ube 065 CordgNj1S) 4 I 


wees Moe GNAIS=dWwilsS OS 


eo! 








a ee 


a ae 


OC 


62€ 


SGe 


Lee 


GGe 


GCE 


HWE 


CCE 


rac 


Oct 


6Tt 


Bie 


ite 


oe 


GTe 


7Le 


lee 


GN3 


qos oo6 


(1° GS] AG! xS "77 Ga = eee G87 co LSS 

S= ONO INIGHG* XS* 77 8c ey ne 8 XS 77 7 Ae Gels Va lade © X9°7 7 77/7 
SIGHC8XO8 S777 7/7 © ° MHC © Gaal me i ete ON Canteen ree ed Ve UO 
AQNNOS Vad HOUVAS Wsy Da dec ec ere Jee ot eaeiiene Nl Cumcumeel son ic ° 91 C 


? 


‘ IV WOLVG LSVDYOS *THOZ*X6*S//SGSSS3NLSIC *VHETSX9S//*§SSVIDNNHLS X9T 


JF SIGHS SX9°7/SOLHCEX9 S&S 7/7 Vad SOD Vie ee ee ee IVWHOd O16 


AS XS ONO IWLOT . 


J 


S1LV IWLGSONO IMS ONO 7a ee ale ely ING SSC Geer yen iC ASE INI dd 606 


666 OL 09 

(T*S4S= AHESXGS//ST°GHS= XHESXG6 

©//*©2*L46= ONOIWILGH6SXS*//*7°2 462 LV IWLOH8*XS*///6AG G3SVATSYHITS 
SXOSS/S// 16 AGHZSX9SS/S// 716 *GHZEX6S//6°GNIM ONSZ HLIMHST&XCTS s/o Le 
ISHNSYHVS SARNDAG HOTS T*24*SY AWVS GNIM 3DV4NNS 4O 1LDAS4R °VHEES X69 
6//©*SYNOH HL¢Z*L4*HLV3G TILNA AWIL *SYNOH HGZ*Z°L4*NOTLSAVHXSHOTS 
‘XZI//S11LNN AWIL *UaLVM NI NVW 4O ZH6Z 916 WOYNS NOILIGNOD YOAIY 


Nase ® CHL CS X6/7 © MeO? yeas Meee NC Colds Ad Ot 


114 





APPENDIX © 
Mrareciuue Uri fg: 


This Appendix is an extract from the study, "Computer 
Program of USCG Search and Rescue Procedures," which was an 
unpublished undergraduate study performed by Cadets First 
Class Robert L. Armacost and Norman T.. Saunders at the United 


States Coast Guard Academy, New London, Connecticut, in May 


1964. 








RLA PARACHUTE DRIFT 
10 Feb., 1964 

After examining the Parachute Drift Table, Figure 6-3 
in CG-~308, it seemed that most of the tabular values were 
linearly related. I plotted these values as a family of 
curves of wind speed. I used the height of the parachute 
Opening as the ordinate, and the actual drift as the abscissa. 
The curves didn't prove to be exactly linear but they were 
very nearly so. For some reason, I don't feel that this 
Linearity should be so apparent. JI decided to try to arrive 
at some purely analyticaliequation to describe the motion of 
the man as he comes down. 

If we assume that he is jumping from a great height, we 
know that he has to free fall to at least 18,000 feet. At 
higher altitudes, ne may suffer from anoxia, or stands the 
chance of freezing since the Petipa soe eee ~5.3°F. In order 
to consider his motion in free Hanae we must examine the com- 
ponents of his motion. They are due to: 

1. The force of gravity. 

2. The effect of the wind on him in a horizontal - 
earection. 

3. His motion due to the velocity of the airplane. 

All of these components are opposed by the air resist- 
ance. The need for these equations is so we may have a very 
accurate entry for the computer. | 


The derivation follows. 
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Vertical Component of Free tall 


Let the y-axis be the vertical axis and be positive down. 


Air Resistance, R = cy. 


Definitions: 


We may 


m = mass of the man 

g = acceleration due to gravity (assumed constant) 
= 32.2 ft/sec” 

W = weight of the man = mg 


Qn = dvm = acceleration of the man 
dt 


Rs air resistance 
2 


@* = factor of proportionality 


now start to write our equations. 


navn = mg - evn oa 
8, 


Let k* = ck, then _1.dvm = o ~ Vin“ 
1 k* dt k 


Let u2 = 2/k* 


dvm = k’dt 
CL =< Vier > 


f 


Vo = O at time t = 0 


t 
[ae rt kat 
_l.tanh— Om ye K*t - got 
[ee Pa 


u 
Vy = uetanh (G-t) (1) 


Py 





AS t>°o, vy, =u as a limiting velocity. 


We may write further equations 
a = uetanh(6°t) 
aw 


a rire, 
oe WF i 0 


At t = 0, 


ig. 

g | 
ay and uetanh( iS et) 
Guts u 

fe) © 
y = _u_In cosh(—B—t) 
g/u ee 
¥ = _v“in cosh(—S—et) 
& 
Let v= _£. 
Uw 


Then: y s _uiIn cosh(\t) 
A 


Av = In cosh(At) 
U 


(A) 


e ~ cosh(,t) 


G = ag soon § 
DN ; 


Refer to page 51, Betz, Burnham, Ewing. 


for a man in free fall. 


mate value of c. 


of motion: 


(2) 


(3) 


They say u = 175ft. 
SEC. 


From this we may evaluate an approxi- 


Se ee 
Jn 








uss/en = VW 
C e 
c= Vy = 4/180 
u 15 
Cua ©0756 (A) 


Horizontal Component of Motion Due to Wind 


Assume: 1. Horizontal, steady wind 
2. The average man is 6 feet tall and wéighs 180 
pounds. He presents a projected area, A, = 1125 in? 
Neglecting air resistance, we can write a few funda- 


mental equations and find an equation for this motion. 


F = force (pounds) 
p = pressure (pounds per square inch) 
An = projected area of the man (ft*) 


Vy, = horizontal velocity of the man 
Viz = velocity of the wind 
F = PAm = Map = mdVy, 
at 
Pz X Vir where velocity head = Ve 
Ze 2E 
and density (pounds /ft?) is equal to ¥ 


My AVs _ alec hee 


at ZE 
But ny = Win 
: & 
Z 
GVn = = Vw Am = _¥ Vy An 
dt ZEN Z\im 


Pig 





be 2 
Let K2 = aon An 
Z Ym 


Then: 


Ke VypofAn_ = V gy REISE 
Vem V2(180) Clad) 


Ke (4.16 x 107-)¥,, (B) 
K* = 17.33 x 107*v,* (C) 
nooks 

ae 


Vn = K*t = 17233 x 1074V;,,*¢ 


For a wind velocity of 60 knots (101 ft/sec.): 


For: 


Vin = (17633 x 1074) (101) 
ce 615 £t/sec., R= cv 


15 ft/sec<V < Sonic, Re eX ye 


To get some idea of our proportionality, Iet us assume a 


free fall of 2000 feet; we must compute the time and thus, V,. 


y = 2000 | 
c = 0.0766 c* = 58.8 x 1074 
m = tn = — = 5.59 

= c2/m = 10.5 x 107% 

u- = g/k* = 3.06 x 10% 

» = g/v = 32.2/175 = 0.184 


Substituting these values back into Equation (3) and solving: 


and 


t = 15,2 seconds 
ve - 17.76 = 269 ft/sec. = 3159.5 knots. 
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Therefore, we may assume that the resistance to his 
horizontal motion is proportional to the square of the 
velocity. Since he is in the same medium, the constant of 
proportionality should be the same since he should have the 
same limiting velocity if sufficient force were applied. 


aVin = K7~c2 Vy" 


dt. 
gies _= dt 
Kr 6°V, 
avi, 5 = oat 
-- - Vy, | 
c 


Then: dV, = edt 
B 2 


- Vb _— 
Att =0, Vn = O 


¥ 
iJ ave = (eat 
one “Vn* O 


= | 2 2 
% tanh ~(V,/B) = ct =_Két : 
B : Re 


Vy = Betanh(K*t/3) (4) 


As t increases, Vp, approaches 8B as a limiting velocity. 


B = KK — ae x dona ™ 
Cc We COT 





— 
| ox = Btanh(K*t/B) 
ite 


rZE 





At Cope Ol, pac, (0) 
x = B“in cosh(K*t/B) 


“Ke 
ice 
X= —1,tn cosh (K*t/B) dae MES cosh / i ee Loh 
c cK a5 SNe. 
x = LAn cosh(3.19 x 1077V,,t) (5) 


Cc 
CG-308 gives the values for parachute drift for a wind 


of given velocity, but this appears to be based on the 
assumption that his initial velocity is zero in both the 
vertical and horizontal planes. However, if we are to deter- 
mine the initial navigational error of the aircraft at the 
instant the pilot leaves the plane, the pilot is traveling 
with a velocity equal to that of the perl 
To write an equation for this component of his motion, 

we will consider that there is no wind and the plane is 
flying horizontally at constant speed. 

We will define the y-axis as the vertical and the - 
x-axis as the horizontal in the direction of his motion. 

Since we know that his initial velocity is that of 
the airplane, and there is no force acting on him except 
that of air resistance, he will decelerate to zero velocity. 
And since his initial velocity is so great, we can readily 
assume that the resistance is proportional to the square 


of the velocity. 
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Horizontal Motion Due to Initial Velocity 


k~ = constant or factor of proportionality 
Vy = velocity due to that of the airplane 


Vo = initial velocity 


2 
Lez t + 
(+ 4.) 
is Vu fob + i 
Vom 


Vok*t + m (6) 
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a 


F 
m In (Vok*t + al 
. O 


aS 


a 
t 


min(Vok*t +m) - me) 


eee + 4 


n Vok* + 
m 


ras 

qt 
f “ns * 
nN 

be fs 


rs 
it 
ct 





1“ 


Le 


Let b 


m 
_— 


7 ia be : 1 (7) 
va 


a 


Here 1 must state that a few other assumptions have 
been made. These are: 

he 6Lhe rotation of therearth Hee USED neglected. 

2. The coriolis force has been neglected. 

3. The sohereicity of the earth has been neglected. 

I feel that these assumptions would not materially 
affect the accuracy of the drift. This is So because of the 
short distances involved. 

Initially, the formulae were derived for conditions 
giving a limiting velocity of 175 ft/sec. This value was 
taken from Differential Equations with Engineering Applica-~ 
tions, page 51. We know that air resistance or parasite 


drag has the form D = Caf A (Theory of Fljght--von Mises-- 


fe 


page 95). This corresponds to c” = tgif A in) OUn ini Glau 





derivation. There we assumed that the value of c* was con- 
stant. A quick glance at this equation immediately reveals 
that c* is not. constant Since it depends on the air density 
which varies with altitude. 

The area, A, is the projected or frontal area of the 
body under consideration. In order to get the maximum 
value of drift, I have assumed that the projected area of 
the man will be the effective area in all three cases of 
motion. Since a body will move with the frontal area normal 
to the direction of movement, we will use this in free fall, 
motion due to the wind, and motion due to initial velocity 
and thus have the maximum drift. This value will then be 
a constant. 7. 

The density of the air, P, depends on the barometric 
pressure and the temperature of the air. Since each of these 
vary with altitude, the density in turn varies with altitude. 
I will use the values of pressure and temperature for the 
U. S. Standard Atmosphere and plot these values against - 
altitude. Hopefully, I wiil be able to write an equation 
defining this variation. 

The third factor is the drag equation is Og, the co- 
efficient of drag. Cg depends on the shape of the body and 
the Reynolds number of the flow. | 

| IR = 2 = Reynolds Number 


a 


=~ velocity 
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tC 
i 


= length 


/._ x kinematic viscosity of the air 


~ dynamic viscosity of the air 


fi 


‘sy 


density of the air 


fr 


The dynamic viscosity, uw, varies with the temperature 


but not with pressure. (Applied Fluid Mechanics - page 338). 
Therefore we may plot “# against temperature and in turn con- 
vert it to altitude. With this and the plot of density, we 
may plot the kinematic viscosity, 7, versus altitude. 

With this, we may investigate the variation in TR 
and thus the variation in Cg. 

Data for plot of Dynamic Viscosity versus Temperature. 
Ref: Table 52 - Handbook of Meteorology.- Berry, Bollay, 
and Beers ~- McGraw Hill ~- New York - 1945. 


Temperature Temperature Dynamic Viscosity Dynamic Viscosity. 


OCG oy eran/ gaps pound~sec,/Pts 
-70 “9k, ee a2I9D : 
50 ~58 147 BOT > 
~30 —22 152 Beso 
-10 th OZ 3.380 

0 32 172 3.595 

10 50 ae 3.700 

20 68 181 Ba vitsiO 
30 86 186 3.880 

LO LOL, 191 3099 
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Data for plot of Barometric Pressure, Temperature, and 
Dynamic Viscosity versus Altitude. Ref: Pressure and Temp- 
erature - Handbook of Meteorology, Dynamic Viscosity - Plot 


of vs. Temperature. 


Altitude Barometric Pressure Temperature Dynamic Viscosity 


Ft.x 1000 mm. of Hg OF pound-sec,/ft* 
(x 1077) 

0 760.0 59.0 3.735 
1 732.9 55 0h 3.720 
2 706.6 51.8 | OH 
3 681.1 LS 4 3.695 
L, OpToPes! Lie8 3.680 
5 632.3 ee 3.663 
6 609.0 37.6 36643 
7 586.4 34.0 3.618 
56k 30.6 3.579 
9 54362 27.0 3.526 

10 522.6 2304 3.473 ae 
il 502.6 19.8 B80 
12 183.3 16.2 3.396 
13 Ob 5 12.6 3.369 
lb Ll, 6 oh 961 3-348 
pis 428.8 565 3.325 
16 411.8 WS) 3.304 


39503 ~1.7 3286 


i 
~] 








Data continued: 
18 
im 
9) 
Cast 
22 
23 
2h 
PAS 
26 
27 

28 
29 
30 
eee 
32 
33 
34, 
Be 
36 
Bs 
38 
39 
40 
Ly 


379.4 
364.0 
349 01 
Bo ae 
320.8 
307.4 
29h oh 
282.9 
269.8 
25S sm 
24.6.9 
236,40 
22586 
Zed 
205.8 
196.4 
187.4 
178.7 
170.4 
N62eu 
154.9 
Zee 
Lhe, 7 
Thee 


ene 

-8,7 
~12.3 
~15.9 
ete) 
~23.1 
~26.5 
-30.1 


= 3.3 67 7 


“373 
~L0.9 
Lhd 
“h7 9 
Pere, 
D5 el 
=O ef 
~62.3 
~65.7 
~67.0 
~67 20 
-67.0 
-67.0 
-67.0 
-67.0 


30226 
3-248 
3.4230 
Beei2 
3193 
32174 
Delle 
36140 
Belee 
Dalle 
3100 
3 4090 
3-083 
3.0756 
320750 
3.076 
BeOSn 
3093 
3-101 
Seth ep! 
3 e1Gr 
3eLOl 
3101 
34101 





Data continued: 


1,2 127.9 
1,3 122.0 
LL, lpley,s 
LD LLOgs 
L.6 105.7 
L7 100.7 
Le 96,05 
LO 91.57 
50 87.30 


~67.0 
-67.0 
-67.0 
-67,.0 
-67,0 
-67.0 
-67.0 
~67.0 
-67.0 


Broken! 
32101 
32101 
Pre LO 
3.101 
3.101 
peel 
3eLO1L 
3 101 


The density of the air and the kinematic viscosity 


may now be calculated. 


p> abe 


where: p = pressure 
acceleration due to grav- 


& 
ity (assumed constant) 


R «= gas constant for Air 


temperature 


aw 

—_ 

ry 4 . > ' 
' ® ¢ 


‘yh 


/ 


Altitude Density = iil Kinematic Viscosity 
Ft. x 1000 Slugs/fs?2 FtS/sec. x 1074 
Ga 10-3) 4) 

0 2 pXe) ILS ASS, 

1 2.005 1.614 

Z ee el,5 eooe 

3 22159 Le yallas 

L Cape ial is. 1.740 
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Data continued: 


Oo Oo N D 


LO 
ie 
ia 
aR 
Lh 
15 
16. 
17 
18 
iN) 
Zo) 
fale 
ap 
23 
Zu, 
25 
26 
27 
28 


2050 
1.985 
1.930 
1.868 
HOLS 
Ue (ee. 
1.705 
1.648 
is 92 
L548 
1.492 
LALO 
1.400 
1.357 
£309 
aac) 
1,225 
force 
Lelhe 
ied. 03 
1.069 
1,028 
0.994 
0.960 


1.789 
Isles 
iyi? 
1.916 
1.942 
1.976 
ZeOL® 
2,008 
eo Lee 
Zeon 
D0 PAPE 
72 aPadsy 
2o3h8 
ey EAL 
20465 


26995 


2.620 
2.704 
Toys 
2.860 
£e9k7 
3.043 
3130 
pee Y: 





Data continued: 


29 0.921 30354 
30 0.888 3.475 
Bal 0.864, S640 
Bie 525 Bb. 730 
33 0.794 Bo7 5 
34 0.765 L080 
35 GUESS da sFadif@) 
36 0.705 4.405 
oy oO, O07: ig Ock 
38 0.639 L860 
By ©6016 22094 
1,0 0.581 bese 
bi 0.555 - 5.595 
L2 0.529. 5876 
43 0.504 eles 
Lh OTE em 6.450 
15 0.457 6.800 
1.6 0.436 feeed 
17 0,416 72470 
L§ 0.397 POSSE 
L9 0.379 See Oe 
50 2301 §.610 


Now that we have the various graphs plotted, we have 


to investigate the variation in the Reynolds Number to deter- 
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mine the variation in the drag coefficient. 


ik = VL 
| V 
At 50,000 feet: = 8.61 x 107% 
= O20) pomhOas 
= Se OL. clone 


At 10,000 feet: 


1.976 x 1074 
22759 x 1077 
Siolayss) Se Shel 


V 
f? 
_ 
VV 
hg 


Va 


The welght of the man need be revised to account for 
the weight of the parachute and the flight suit. If we 
figure 20 pounds each for the last two items, mg = W = 220 
pounds. 

We must assume a shape for the man. There are two 
possibilities: 

we Circular cylinders, @ 2.1.5 feet 

2. Ellipsoid, L-=- 6 feet, d=1.5 feet 


In both cases, the Reynolds Number will be the same. 


At 50,000 feet: Sy en 1.738 x 10°V 
IK 4) 801 arog 


Va 
2, tao LO” 


When the man is wearing a flight suit with boots and 





At 10,000 feet: 7 = Vd 7.59 x 10°V 


helmet and has a parachute on, he will more closely approxi- 
mate a circular cyclinder. For this reason, and the fact 


‘that there is very little data on ellipsoids, I will use the 
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form of a circular cylinder. 

From Weisberger's experiments, if we can keep 
K>5 x 10?, we will have a constant value for ©,. (Applied 
Hydro-and-Aero-Mechanics, Prandtl and Tietjens). 

We will examine JK at both 50,000 and 10,000 feet to 
determine what velocity we must have, 

At 50,000 feet: TR = 5 x 10? = 1.738 x 103V 
V must equal 288 ft./sec. 

At 10,000 feet: ]Z = 5x10? = 7.59 x 102V 
V must equal 65.9 ft./sec. 

To see if this is possible, we must check the limit- 
ing velocity at each altitude. We know that the limiting 
velocity is reached when the drag equals the weight since 
there is no acceleration. - 

With [KR = 5 x 102, Cg = 0.33 for an infinite cylinder 
(Handbook of Engineering Fundamentals - Eshbach). 

Dye we aes 


At 50,000 feet: V* = 2W = 2(220)(1 
CgAA  0,33(0.361 x 107) (11-25) 


v* = 47.3 x 104 
V 685 ft./sec. 





We found that we must have a velocity of 288 ft./sec. so our 


assumption is valid at 50,000 feet. 


he 1OCOO fects av = 2mm 2(220) (uk) 
| | Cg AA ~ 0.33(1.759 x 10-4) (1125) 


T = 9,8 «10° 
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¥ = 314 ft./sec. 

We must have a velocity of 65.9 ft./sec., so our 
assumption is also valid at 10,000 feet. 

This criterion of [> 5 x 10? will hold true every 
where between 10,000 and 50,000 feet. This becomes obvious 
when examining the curves. 

With IR>5 x 10°, we may enter Table III in Rouse, 
page 249, and find that for a length-diameter ratio: of 5 
(which is very nearly what we have), Cg = 0.35. 

Looking toward the same problem with a parachute, 
this investigation has also revealed a value of Cg = 1.33 
for a parachute. (Theory of Flight - von Mises — page 98). 
This value is in agreement with the values found in other 
texts for a hemisphere, hollow eR eneatee The criterion here 
is that TR > 102 which is easily recognized. 

It should be noted that these values of critical 
Reynolds Numbers decrease with increasing turbulence of the 
air. The values used herein are for air with no turbulence 
in the approaching stream. (Fluid Mechanics - Cox and 
Germano ~- p.242). This gives us a factor of safety since 
the air will be turbulent. 

We have now covered the factors with one exception, 
density. Therefore, the problem at hand is writing an 
equation which describes the density variation. 


From the plot of density versus altitude, it appears 
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that we have a decreasing exponential. This will give us an 
equation in the form of 
fo = heF 

A plot of vs. altitude on semi-log paper will give 
us an indication of whether or not we have an exponential 
relationship. (See graph next page). 

By sighting down the plot, it is seen that we do not 
have a straight line. However, a straight line fits the 
curve from O to 32,000 feet, and another from 32,000 to 
50,000 feet. 

if would rather avproximate the relationship with one 
smooth curve rather than make it discontinuous. The problem 
we now have is that of tolerances. How much tolerance can 
we allow and where should errors be a Roscmne 

With this determined we only have to integrate and 


moe man's motion is Gefined. 
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RLA 
11 March 1964 
From the plot of VS. y on semi-log paper, we may 
approximate the curve with a straight line. ‘The equation 


has the form 


let = 2.40 x 107? 
For our first approximation, consider 

y = 35,000 feet, where z= 0.735 x 1073 
From this, B = -3,375 x 107) . 

With this as a starting point, we may find the dif- 
ferences between the formula value and the actual value. 
The value of B may then be adjusted to give differences which 
are acceptable. 

Several adjustments give B = -3.44 x 1079. For this 
value, the following differences exist. 

B= -3.44 x 107? 

(Values for Actual density, formula density and the difference 
between formula density and actual density are all exvressed 
as slugs/ft.? x 1077. Altitude is in thousands of feet). 


Altitude Actual Density Formula Density Difference ° 


O 2.38 2.40 0.02 

5 2.056 2.021 =O. O20 
10 1.759 3 1.700 -0.059 
15 1.492 1.433 -0.059 
20 : 1.264 1.209 0.055 
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Data continued: 


25 1.069 1.040 -0.029 
30 0.888 0.856 -0.032 
35 0.735 0.731 -0.004 
£0 0.561 0.607 0.026 
45 O457 Oe ALC - 0.053 
50 Oecoll 0.430 0.069 


We now have an equation for the density of air as a 


Mmunctsion Of altitude. 


(P= 2.40 x 1073e3 oh x 1O™y) 

Now we can proceed and integrate our basic equation 
for free fall. At this time, I will use the letter h to 
designate altitude and y will be used strictly for vertical 
displacement. Therefore, y = ho - h where ho is the alti- 
tude at which the fall began. the relationship for air den- 
Sity now becomes | 

AP = 2.h0 x 1073e!-3 hh % 107h) 


Gur primary equation is now 
=> maV = mg - WA A 
dt 2 
let CqgfA, = C-. Then mdUp= mg - CVZe™ 
Z at 
Since ysh)o -h, hz=hg - y where y = {Vat 


ndvV, =x Mg - C ie e Bho 
re B 
‘eh eos 
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Bho . 
° is constant. 


For a given altitude, e 
Let C.Pho ome, 10 
Since 8 is negative, we have 


mV = mg - DVZe"Y 
at. 


Ll have tried for over a week to integrate this equa- 
tion consulting every Aipcorentsad equation Leb Oe Chat 
could find. This has not vroduced any positive results. 
Since this equation appears to be impossible to integrate, 
we can also approximate the relationship for density as a 
series of straight lines of changing slope. We may consider 
each curve over a definite interval and take three intervals 
in order to make a good approximation. This equation has 
the form as. 

Pp - an +b 
The first interval we shall consider is 

GQ) oe ws 14,000 feet 

feo 2 Gwe) 2057 oe UO 

7 by, = 2:37 =e 

At h = 14,000, © = 1.53 x 1073 

=> a; = -6.00 x 107§ 

Therefore, Oz = -6.00 x 107%h + 2.37 x 1079 


The second interval is 


IAN 


30,000 feet 
At. h = 1h,000, 0 = 1,53 x 1073 


(2) 14,000 <b 
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=> boa Weapon 

At h = 30,000, = 0.87 x 1073 

ede ee ome 10-8 

Teena P2 > whel25 x 107°h + 1.53 x 107 
The third interval is 

(3) 30,000 $n 

ft h = 30,000, P = 0.87 x 107? 

=> b3 = 0.87 x 107? 

At h = 50,000, PP = 0.34 x 1073 

=> a3 = -2.65 x 1078 

Therefore, [? 3 = =2,05 2 1078h + 0.87 x 1073 


iA 


50,000 feet 


With our density described in another form, we may 
again write the differential equation for the man*s motion. 


Now we have 
| fe 
mdVy = mg - Cglah + b)VEA 
= 5 1 
mdVn = mg - CgakehVG - CabinVi 
Ze Zz 
Let Coin = C,, then mdVm = mg - C.,ahVA - C bY Es 
a dt 
But again, ¥ = [Vdt =hg- h, andh =hy - y 
mdVm = mg - Caalhg - y)Vq ~ C_bVii 


Tsaes 
mdVn = mg - (C,ah, + C,b)VA 4+ CayVS 

eo & aw"O a nt} qty *m 

Again, I have worked several days on this equation and 


have been unable to integrate it. 


if we can’t integrate this equation, or the one 
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involving the exvonential relationship for joubity, then we 
must find some way of approximating this solution. 

Assume for the present that we have a constant density. 
In that case, our original integration holds and we have 

Cee bane A 

(2) y = Cu/A}In cosh(n tt) 

2 ee = Lcosh™ dv/t ) 


ti 


ale Ns Vf Bes 
Bog PAy aw 


Therefore, u = 2 


In order for us to test the validity of this approxi- 
mation, we must determine the numerical values and express 


these constants as functions of the density. 


u=_/_20 = 12.7/Vp 
Cg P Fan 


= 2535/0 


ab 


14] 





ew Ss 5402 
fo alr 2W = 0,394 
& Ca PAn {PP 
ceil in = 0.199/ 


i 


|» aa >: 


We may now rewrite our equations as 


4 (1) Vin = 12sZtanh(2.535-/ %) 
~/ 
(2) y = 5.02 1n cosh(2.535-V 0 t) 
-/ P 


(3) te 0.39nc0sh7}|¢(0-199,0 y)| 
~{P- 


If we analyze this method of solution over a large 
vertical displacement, we find that effectively, we are con- 
Sidering the density to be constant in accordance with our 
initial assumption. This will provide very inaccurate 
results. Therefore, if we are to assume constant density, 
we must limit our value of y to a small increment. I will 
arbitrarily choose a value of 1,000 feet for y as the inter- 
val of disvlacement. This I feel is a safe assumption and 
should not introduce any significant error. If we are to 
proceed in this manner, we must account for the initial 
velocity at the beginning of each of the intervals. To do 
this we must return to our original derivation where 


miVn = mg - ¢*Va 
dt 
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BAS = ce Vin 
ho t = O, V, 


Tih. 


=. Vp 


If we integrate, 


. Ng 
tanh“2(T/u)} ~- (g/u}t 


0 


tanh1(V_/u) ~ (g/u)t + tanh71(Vo/u) and 


Vy = utanh|(g/u)t + tanh (Vo/u)| 


Integrating again 


Y= 

6) ga 
mn 

t = 0,39 
ay 2 


-(u*/e)In cosh{(g/u)t + tanh-1(V,/u)] 
[cosh | ~ tanh™(Vo/u)| 


oen"2( 02-2957) ~ vant” Y(Yaxl@ 
12 7 


To pursue this method of solution, we must use the 


following metnod of approach, and we may use the exact 


values of air density. We will consider hg = 20,000 feet. 


For y = 1,000 feet, we will determine the time for free fall 


to 19,000 feet.assuming a constant air density equal to the 


value for 19,000 feet. We will also determine the velocity 


at 19,000 feet which will be the initial velocity which we 


shall use for the interval between 19,000 and 18,000 feet. 


Results of these calculations follow. 


h 
19,000 
18,000 
17,000 
16,000 


px 107? 
1309 


£6357 
1.400 


LALO 
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VEC: 


Zak 
Zee 
Bie: 


Zee 
fraps 
Jib fe ge 


SeL6 
4 okS 
he 5O 
bedi 





In this case, t does not give the time ie free fall 
except for tne first 1,000 feet. It appears that here the 
man has reached the limiting velocity. Since the velocity 
changes very little, we may divide the average velocity for 
the interval into 1,000 feet to find the time of fall for 
that interval. 

If we started the same procedure at say ho = 25,000 7 
feet, and at 17,000 feet had a velocity of 222 ft./sec., we 


could be reasonably sure that our approximation is correct. 


Let hy = 25,000 feet. 


: px 107? Va Vy t 
21,000 1.103 227 0 Saue 
23,000 1-142 228 227 ees 
22,000 1.180 226.5 228 Lom 


The tendency of these calculations is to verify each 
other. Our first thought would then be that this approxi- 
mation is correct. But, let us look at the limiting veloci- 


ties at the different altitudes. . 


h - 17,000 feet 
v .  2W 
Cg P Ay 
V = 339 ft./sec. 
hs 2h, ,000 feet, V = 382 ft./sec. 


By our equation in the approximation, consider the 


‘interval of fall from 50,000 to 49,000 feet. 





49,000 £h $ 50,000; = V = 245 ft./sec. 

48,000 $n £49,000; V¥ = 228 £t./sece 
The limiting velocity at 48,000 feet is: 

V = 636 ft./sec. 

In all three cases, we notice quite a disparity be- 
tween the velocity as determined by our equations and the 
Pimiting Veloci tweens irom bois, = & think we must accept the 
approximations as being invalid. 

In order to determine the time of free fall, another 
method of approximation could be used. We could assume that 
the man has reached his limiting velocity and plot these 
values of limiting velocity versus altitude. By taking the 
average velocity over a 1,000 foot interval, we could find 
the time taken to traverse that ere Now, suppose that 
this would give us results without appreciable error. We 
would then have no means to determine his horizontal motion 
because we have no knowledge of whether or not he will reach 
his limiting velocity in the horizontal direction. : 

For an accurate Sethian then, we must integrate our 
equations of motion. JI will now write these equations and 
define all of the terms in summary form. 

Our first set of equations, (1), (2), and (3), are for 
the condition of free fall. Bear in mind that the equation 
for the density of air may be either the exponential rela- 


tionship or the linear relationship. 
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(1) Motion vertically downward. 


e 
maVm = me ~ Gq {foo Vien 
at 2 {ah +b 


. (2) Horizontal motion due to wind. 


oe nn 
2 
maV, - a A aes VE 
dt - - og ah +b am 
(3) ee motion due to initial velocity. 


gu - - _ 2 Linn 
an +b 


It is easily seen that the results of equation (1) 
relating the displacement in free fall must be an input into 
equations (2) and (3) in the relationship for density. 


The various terms are defined below. 


m - mass of the man = _W_ in slugs 
& 
W - weight of the man in pounds 
& = acceleration due to gravity = 32.2 ft ./sece” 
Vm = velocity of the man in the vertical direction in 


—_ 


ft./sece 

Vi = velocity of the man in the horizontal direction 
due to the wind in ft./sec. 

Vi, = velocity of the man in the horizontal direction 
due to the initial velocity of the aircraft in 
ft./sec. 


Cq = coefficient of drag = 0.35 
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vrojected area of the man = 1125 square inches 


Vw 2 velocity of the wind in ft./sec. Bh 
e 
\ = Welehe dens Ua oreume sii =p f= 5 Pee 
| ah + b 
Po = PAPO) sic 1077 slugs/ft2 


B= -3.44 x 107? fttt 
nN = altitude in feet 
Ao = altitude of bail-out 


y ~ vertical displacement in feet 


a and b are defined for the following intervals: 
o $n $ 34,000 a = -6,00 x 107° 
b - 2 oa x on 
Te 000 < ty S 305000 MINE ay =eenti2> co 10m 
b= 1,53-x 107 
30,000 $h £50,000 =a = -2.65 x 1078 
b 2 @687 x 1072 


We must now Look at the equations for fall with a 
parachute. Basically, there is no difference except for the 
values of ©, and the projected area. : 

Mons ROO ua the vertical direction, we may use a 
value of Cg = 1.33 (Theory of Flight - von Mises - page 98) 
for a parachute. This value is in agreement with those | 
found in other textbooks for a hemisphere, hollow upstream, 
for a Reynolds number greater chen 107, Hereafter, I will 
refer to this drag coefficient with the symbol Cays 


Since we have assumed the shape of a hemisphere, the 
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projected area for motion in the vertical direction will be 
the area of a circle of a radius equal to fourteen feet. 
Therefore, Ay = (14)? = 615.u ft? 

For motion in the horizontal direction, the projected 
area, An = bath)? = 307.7 ft.” This is the area of a semi- 
circle of a radius of fourteen feet, 

The determination of the drag coefficient for side- 
wise motion presents a problem Since there have been no 
experiments performed for a hemisphere with its major axis 
perpendicular to the direction of flow.: In the vertical 
direction, Cgy = 1.33 with the hollow end upstream. For a 
hemisphere with its hollow end downstream, ©g = 0.34. After 
checking in many fluid mechanics textbooks and studying the 
various drag coefficients for different bodies under several 
conditions of flow, and studying theoretical flow patterns 
and pictures of actual flow patterns, 1 will assume a drag 
eoefficient for the parachute in the horizontal direction. 
The value which I have chosen is Og, = 0.60. . 

Nore we may write our equations of motion for para- 
chute fall. 


{i) Vertical nia 


mdvm = mg - Cay Vek, 
at an b 


(5) Horizontal motion due to wind. 


€ 
méTn 2 XVSkn - CapdMPor Cv2A, 
dc. 2g ah + b 





(6) Horizontal motion due to initial velocity. 


bern 
C 2 
maVy = -’dh or b Ay, 


dt 2 (ah +b} ~ 
tie Ce 
Ay = 615.4 fs 
Cah = 0-80 


Mn oO) Stes 

The rest of the terms are the same as defined for 
equations (1), (2), and (3). 

In all of the above equations, all must be integrated 
twice in order to give a displacement in the particular 
direction. Equations (1) and (4) must be solved for t by 
first integrating to find the displacement as a TunclLon om 
time and tnen transposing to find time es a function of dis-~ 
placement. i1f we call the time of free fall given by equa- 
tion (1) ty, then t; substituted into equations (2) and (3) 
gives the displacements in the horizontal direction. In the 
integration of equations (5) and (6), consideration must -be 
given to the fact that there is an initial velocity. This 
initial velocity is the value given by the first integration 
(for V) of equations (2) and (3) at time t,. 

By summing up all of the horizontal values considering 
the direction as well as the AVsikkande given by our equations, 
over the period of time from tg to ty (free fall) and from 


ty to ty (parachute tale we have the parachute drift. 
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L would now like to write this problem in the form of 
&@ program flow chart. When the six mits of motion are 
integrated, these may be written as additional sub-routines 
and added to this program with a predefined connector. I 
will also omit all machine instructions since all processes 
are mathematical operations. 

A few new terms were added in the flow chart which 
snould be defined. 

A= the direction in degrees toward which the wind is 

blowing. 


the course of the aircraft when the man bails outs. 


Ul 


p 
O 
After all of this work, I feel that I should explain 


v{ 


“whe resultant direction of parachute drme. 


why I have put this much time in what ie seunia@ir 2 minor 
matter. ina recent study of the feasibility of the use of 
a personnel locator beacon, the following facts were brought 
to light in a study of Naval aircraft accidents. The study 
showed that over 98% of all incidents occurred within 60. 
miles of a rescue facility. This is within the range of a 
good DF fix. Assuming that we know the location of the air~ 
craft when the person leaves his eens, we should be able to 
find him as quickly as possible. By accounting for the 
variables discussed in the analysis of the problem, we may 
determine where he will be in the water with reasonable 


accuracy. in cases where the sea and air temperature are 
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ritical, or where the man is injured, it is absolutely 
necessary that we rescue him as quickly as possible. It is 
apparent that this same solution is applicable over land 
areas. 

This analysis of his motion in the air will make this 


possible upon the integration of his six equations of motion. 
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Chart ID AB ' Chart Neme: Parachute fall 
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Chart ID AC Chert memo: Perechute driftedirection 
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Chart ID AD 


Chart neme: Parachute driftemegnitude 
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APPENDIX D. 


EXTRACTS FROM THE STUDY: 


"PREE FALL AND PARACHUTE DRIFT" 
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APPENDIX D. 


EXTRACTS FROM THE STUDY: 


“FREE FALL AND PARACHUTE DRIFT" 


“Free Fall and Parachute Drift" was an undergraduate study performed 
by Cadet First Class Darvy M. Cohan at the U. S. Coast Guard Academy, 
New London, Connecticut, in the spring of 1965. Cohan expanded and 
refined the work done by Armacost and Saunders which appears in the 
preceding Appendix. Cohan expanded the problem to include the effect 
of winds at three different altitudes blowing with different forces and in 
different directions. He also considered the possibilities of the 
parachute opening at three different altitudes: (1) the altitude at which 
the aircraft was flying, (2) automatically at 14,000 feet, and (3) at the 
surface of the Earth (parachute fails to open). 

Cohan developed three computer programs, one mainline program 
and two subroutines, all of which are included in this Appendix. A brief 
description of each program precedes the program itself. The flow 
diagrams are taken directly from the original paper. 

6 
VECTOR SUM PROGRAM 


Cohan's entire program hinged upon the development of a program 
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to sum vectors in three dimensions. The path of motion represented a 
resultant vector which made an angle 9 with the horizontal. Using sine 
and cosine functions of that angle he determined the vertical and 
horizontal distances a man would travel in a specified time interval, 
given the velocity component, V, within the air mass. The total 
horizontal displacement relative to the air mass is the sum of the 
distances for all time intervals and is represented by DTR. By applying 
the horizontal motion of the air mass the final resultant motion with 
respect to the Earth's surface is obtained. 

The theory behind his vector sum program, VSUM, was this: 

1. Establish one vector as the major axis used by the computer. 

2. Build a triangle using the angle between two vectors. One leg 
would then be the first vector plus the projection of the second onto 
the major axis. The other leg would be the projection of the second 
vector upon the orthogonal axis. Dividing the first leg into the second 
leg produces the tangent of the angle of the resultant vector. 

3. Design a dummy system to work with the azimuth angles found in 
the real world. Convert the azimuth angles to the dummy system, solve 


the problem, and convert back to azimuth angles. 
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SUBROUTINE ORDER 


In considering the three possible parachute-opening altitudes, 
Cohan ended up with three possible vectors at the Earth's surface, DISA, 
DISB, and DISC with their associated angles, ANGA, ANGB, and ANGC. 
His SUBROUTINE ORDER takes the three vectors, arranges them from 


maximum to minimum and then produces the vector average of the three. 


MAINLINE PROGRAM 


The mainline program, called SUBROUTINE PARAD, uses the two 
preceding subroutines and the basic equations found in the study by 
Armacost and Saunders. The Euler method of statistical integration was 
applied for time intervals of two seconds in the free fall phase and four 
seconds in the parachute drift phase. The exponential form of the 


equations involving air density were employed. 
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APPENDIX E 


A Paper: 
COMPUTER PRODUCED SYNOPTIC ANALYSES OF SURFACE 


CURRENTS AND THEIR APPLICATION FOR NAVIGATION 


by 
W. E. Hubert 


Commander, United States Navy 
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ABSTRACT 


The available methods for estimation of wind currents, mass 
transport velocity by waves and permanent flow (thermohaline 
gradient current) are briefly summerized and a simplified computer 
approach is outlined. 

The computed synoptic surface currents are compared with monthly 
mean current charts and with surface wind conditions. This analysis 
indicates that the surface currents are greatly wind-driven. A detailed 
verjfication procedure which will use the observed changes in sea 
surface temperature is outlined. 

The use of the synoptic current fields for computation of divergence 
and convergence and the resulting changes in subsurface thermal 
structure is described. The relative importance of the synoptic surface 


currents in ship routing, rescue operations and other practices in reviewed. 
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1, INTRODUCTION 

A number of naval, fisheries and other martime operations re quire 
a knowledge of the direction and speed of surface currents, as well 
as their past and near future behavior. The Fleet Numerical Weather 
Facility (FNWF) at Monterey, California, became interested in ocean 
currents primarily because of their importance in Anti-Submarine (ASW) 
applications. 

Large variations have been observed in thermocline depth which 
cannot be explained by mechanical or convective mixing. These changes 
exhibit cycles which correspond closely to the evolution of synoptic 
weather patterns over the ocean. It is quite clear that current atlases 
cannot be used to predict thermocline depth when considerable change 
Can occur in a period of a few days. What is needed for this particular 
problem are daily current analyses and prognoses. 

Navigators also undoubtedly find that atlases, monthly mean charts 
and the like frequently do not give an exact enough answer to the question: 
what is the current at a given point in space and time? The purpose of 
this paper is to report on an attempt to compute surface current flow 
on a quasi=synoptic schedule and to show some preliminary resulis. 

If these results appear to be of use to navigation, means will be found 


to accomplish dissemination. 


ae, 








2, BACKGROUND 

Before selecting and justifying an approach which would be simple 
but still give some hope of yielding useful results, it was necessary 
to screen and evaluate a voluminous amount of literature on the subject. 
Fortunately, a great part of this review had been done recently by 
Laevastu (1962). 

Most theoretical approaches have been mainly concerned with 
explaining the general, more permanent features of the horizontal 
circulation patterns (see, e.g., Robinson 1963). The Russians have 
recently applied correlation theory in an attempt to forecast detailed 
current changes from a known field; however, our Knowledge of the 
initial state (and particularly its derivatives) is often rather sketchy. 
Actually, the ocean responds quite rapidly to hourly and daily changes 
in driving forces, and currents are known to be variable in space and 
time (Knauss 1960). 

These considerations dictated use of a method which would account 
for fairly rapid response and would stand up to daily verification. Many 
aitempts at Current prediction have been disappointments because of 
overSimplifications resulting from the assumptions made. It seems 
logical to separate the total current into its elementary components to 
see which should be neglected, which can be simplified, etc. This is 


the attack which has been followed in this investigation. 
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3. COMPONENTS OF SURFACE CURRENTS 

Surface currents are caused and influenced simulianeously by a 
number of forces which vary independently from each other in space 
andtime. If one neglects the special effects due to variation in 


denth, coastal configuration, runoff, etc., the current vector ata 


given location, time and depth below the surface Deer can be 
given as the resultant of the following components: 


where Weg is the permanent flow (chermo-haline gradient current or 
“characteristic current" as used by Palmen (1930) and Hela (1952)), 
Ww is the current due to transport by wind and waves, W, is the 
periodic part of the inertia current and W, is the periodic part of 
the tidal current. 

The Computations which will be described here cover a period of 
24 hours, and it will be assumed that semidiural and diumal tidal 
components will equal out and can be neglected. In addition, inertial 
eddies will not be considered because the available quantitative infor- 
mation about their behavior does not warrant their inclusion in this 
Simple technique. The two components which this study will attempt 


to evaluate are thus the "characteristic" or permanent transport and 


the transport due to wind and waves. 
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4, THE CHARACTERISTIC COMPONENT 

The characteristic cbmponent is directly related to density 
gradients caused by areal differences in heating-cooling and 
evaporaiion=-precipitation. Although what we usually call “the 
permanent flow" is strongly conirolled by the large-scale, more- 
or-less stationary wind systems, only the thermopmhaline influences 
are included in that component here. Wind and wave effects will 
be lumped into one computation to be discussed later. 

Several workers have found (e.g. Yasui 1957) that there is a 
close correlation between ocean temperature distribution and 
dynamic depth anomalies. Neglpcting salinity, one can apply 
the well-known metéorological thermal wind relationship in the 
ocean if one knows the mean temperature of the layer between the 

Surface and some level of zero current velocity. The characteristic 
current is then given by: 
= 942 oT x K (2) 


Wo ay 


where g is gravitational acceleration, f is the Coriolis parameter, 





T is the mean temperature above the level of zero current, ws 


is the depth to zero current, and [fK is the unit vertical sector. 
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The determination of representative mean temperature (T) 
is, of course, the critical factor in this part of the problem. The 
temperature structure of the ocean is certainly not constant, partic- 
ularly closer to the surface; therefore, semi-synoptic temperature 
fields should be used if possible. The only place where sufficient 
data is available for reliable analysis on a daily basis is at the 
surface so it was decided to use the FNWF Sea Surface Temperature 
(SST) analyses based on 84 hours of reported ship engine injection 
temperatures at the top of the layer. In order to include a part of 
the deeper temperature structure, the SST field is presently combined 
with a climatological field at 209 —_— depth to obtain 
(3) 


T en ok 


Ievate 21200 


Finally, this field is modified empirically in areas where salinity 
considerations are known to be important (Cyashio, Greenland, 

Labrador currents). This in effect corrects the ocean temperatures 

for salinity much as _ the meterorologist corrects atmospheric temperatures 


for moisture content when he uses the concept of "virtual temperature. " 
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5. THE WIND COMPONENT 

According to Ekman (1906), the direction of the wind current 
at the surface is 45° to the right of the wind in the Northern 
Hemisphere and this angle increases with depth. Recent inves- 
tigations reveal that the deflection is more nearly 12-20 degrees, 
being larger and more irregular at.the lower wind speeds (possibly 
because of the increased importance of other components) and 
smaller at higher wind speeds. As the surface wind is about the 
same angle to the left of the geostrophic wind, it is assumed herein 
that the direction of the wind current is the direction of the geostrophic 
wind. 

Numerous empirical studies have indicated use of a single factor 
to relate surface current speed to wind speed. The formula of 
Witting (1909) appears to agree well with available data and further 


allows approximate incorporation of mass transport by waves ina 
simple expression: 


Wy, = Kz 7 W, (4) 


where Wg is the mean geostrophic wind speed for a 24-hour period. 
In the present work it is assumed that the current is relatively 
uniform and unidirectional in the turbulent mixed layer down to the 
thermocline (or about 200 meters). The mass transport of the waves, 
however, would modify this picture as it decreases exponentially 


with depth (Masch 1962). Therefore, if W.. is in meters/sec and 


g 
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Wy an cm/sec, K is taken to be 3.3 for surface currents (ship routing 
and drift computations) and 2.2 for the average current down to the 
thermocline (convergence/divergence computations). Obviously, 
there is a time lag between the change of the wind and response of 
the sea. This lag seems to be shorter than previously believed, 
however, and is partially minimized by the 24-hour averaging. 

Since all computations are carried out in the standard FNWF 
grid system, u and v current components are determined at approx- 
imately 200 nautical mile intervals for all Northern Hemsphere. | 
ocean areas. From these components, direction and total transport 
(nm/day) fields are determined and stored on magnetic tape for later 
output in chart form or as special messages giving the currents at 
specified latitude/longitude intersections. 
6. RESULTS 

Figure 1 is a hand analysis of one of the first current computations 
made on a synoptic schedule (18 GMT 16 November 1964). The con- 
tours represent total current transport in nautical miles per day; 
direction arrows have been plotted in the most significant current 
systems. One can clearly distinguish such well-known features as the 
Gulf Stream, Sargasso Sea, Labrador Current, Kuroshio and Cyashio. 
The low=latitude, westerly return flow which results primarily from 
the “wind component" term is well defined in both the Atlantic and 


Pacific. A narrow equatorial countercurrent was obtained as a result 
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of the 200-meter temperature structure used in the characteristic 
component. 

Figures 2a and b are synoptic current charts for the Pacific and 
Atlantic, respectively, which have been drawn automatically on an 
incremental x-y curveplotter. Each chart requires approximately 
one minute to complete and is of sufficient quality that it can be 
used immediately for radio-facsimile transmission. 

Figure 3 is a climatological current chart for winter. It can be 
seen that many of the most impprtant features are correctly depicted 
in this approach both in location and intensity. 

The problem of automatic plotting of direction arrows on these 
charts has not yet been solved. However, a possible substitute 
has been found and is now being programmed for numerical testing. 
Since u and v current components are available in field form at all 
grid points, it is believed a stream function field can be determined 
iy a relaxation solution of the Poisson equation 
ay ou 
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Vv = - (5) 
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This would permit plotting of a second set of lines (WY) which 
would every where parallel the direction of flow. 
7. VERIFICATION 

A synoptic current chart would be of little value if it could not 


be verified and the computational scheme tuned as required. ODirect 
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current measurements in the open ocean are too few and drift 
calculations made from navigational fixes are frequently inaccu= 
rate in weak current areas, so it is difficult to make a direct 
evaluation. It has been necessary, therefore, to resort to indirect 
means which are susceptible to verification on a synoptic basis. 

Sea surface temperature (SST) is the only oceanographic element 
which permits a reasonable complete synoptic analysis on a hemi- 
spheric scale. Such analyses are made twice daily at FNWF 
Monterey (Wolff 1964), and their resolution is such that SST temper- 
ature changes can be determined for periods of 24, 48 hours, etc. 
From these changes will be subtracted the local changes computed 
from air/sea heat exchange equations. If the remainder correlates 


well with the advective change indicated by W "a VY SST, the 


XYZ 
computed currents can be assumed to be. reasonably correct. 

This method of verification is now being programmed and 
numerical results are not yet available, Subjective study of 
SST change charts and corresponding current charts does, however, 
indicate that the approach described here is useful. It is evident 
that the wind component term predominates in many areas, and 
that it is this term which is mainly responsible for the rapid 
response of sea surface temperature changes in the ocean. 

There are a number of modifications which must undoubtedly 
be made to this program; it is hoped these will be uncovered 


during the verification period. One obvious question is - what 
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effect does thermocline depth itself have upon the surface current 
speeds above the thermocline? 
8. APPLICATIONS 

The surface current program was initiated primarily to determine 
divergence and convergence and the accompanying up and down 
movement of thermocline depth. Furthermore, the results will be 
used for forecasting the advective part of sea surface temperature 
changes. 

Over large parts of the oceans the currents have liitle direct 
effect on navigation. In some areas, however, they should be 
taken into consideration in Optimum Ship Routing. Charts of this 
type should also prove useful in the prediction of ice movement 
and in any rescue operations. 

It is planned to make these Computations on a daily synoptic 
schedule (probably to 06 and 18 GMT), and they could be trans- 
mitted from Fleet Weather Centrals in either facsimile or special 
message format if such is desired. Groups such as the Institute 
of Navigation may determine that there are applications in navigation 


which could make use of these products. 
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APPENDIX F. 


SURFACE DRIFT ELLIPSE PROGRAM 


The Surface Drift Ellipse Program which comprises this Appendix 
was developed by Donald A. Burns of the U. S. Naval Oceanographic 
Office. The specific instructions for using this program may be 
obtained by writing to Mr. Burns, Code 3800, U. S. Naval Oceanographic 
Office, Washington, D. C., 20390. The basic logic of this type of 


program may be found in A Statistical Rose Program, SP=-64, available 





from the Evaluation Branch, Oceanographic Analysis Division, Marine 
Sciences Department, U. S. Naval Oceanographic Office, Washington, 


me. CO. 
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